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Attempts to use coarse-grained molecular theories to calculate corrections to the random-phase 
approximation (RPA) for correlations in polymer mixtures have been plagued by an unwanted sensi- 
tivity to the value of an arbitrary cutoff length, i.e., by an ultraviolet (UV) divergence. We analyze 
the UV divergence of the inverse structure factor S~ 1 (k) predicted by a 'one-loop' approximation 
similar to that used in several previous studies. We consider both miscible homopolymer blends 
and disordered diblock copolymer melts. We show, in both cases, that all UV divergent contri- 
butions can be absorbed into a renormalization of the values of the phenomenological parameters 
of a generalized self-consistent field theory (SCFT). This observation allows the construction of a 
UV convergent theory of corrections to SCFT phenomenology. The UV-divergent one-loop contri- 
bution to S~ (k) are shown to be the sum of: (i) a fc-independent contribution that arises from 
a renormalization of the effective x parameter, (ii) a fc-dependent contribution that arises from 
a renormalization of monomer statistical segment lengths, (iii) a contribution proportional to k 2 
that arises from a square-gradient contribution to the one-loop fluctuation free energy, and (iv) a 
A;-dependent contribution that is inversely proportional to the degree of polymerization, which arises 
from local perturbations in fluid structure near chain ends and near junctions between blocks in 
block copolymers. 



I. INTRODUCTION 

The statistical mechanics of polymer mixtures and 
block copolymer melts exhibit some universal features 
that are well described by self-consistent field theory 
(SCFT). The phase behavior of homopolymer mixtures is 
reasonably well described by Flory-Huggins theory. Var- 
ious inhomogeneous structures formed by flexible poly- 
mers, such as interfaces and ordered phases of block 
copolymers, are accurately described by a SCFT of in- 
homogeneous liquids, which reduces to Flory-Huggins 
theory in the case of a homogeneous mixture. A self- 
consistent field approximation also underlies the so-called 
random phase approximation (RPA) for the structure 
function S(k) in homogeneous mixtures: The RPA is ob- 
tained by using SCFT to calculate the susceptibility of 
a liquid to a hypothetical infinitesimal perturbation, and 
using the correlation-response theorem to relate this lin- 
ear susceptibility to the corresponding correlation func- 
tion. 

SCFT is a highly successful theory, but not a perfect 
one. Among its limitations is the inability of the RPA 
to accurately describe fluctuations very near a critical 
point in a polymer blend or near an order-disorder transi- 
tion (ODT) in a symmetric diblock copolymer melt. The 
range of temperatures over which deviations from the 
RPA are significant is believed to decrease with increas- 
ing degree of polymerization N: The fractional width of 
this so-called Ginzburg region is predicted to decrease 
as A^ 1 with molecular weight JV in a homopolymer 
blencU^i 3 -, and as A -1 / 3 in a symmetric diblock copoly- 
mer meltj^ For molecular weights typical of experiments, 
however, fluctuation effects that are ignored by SCFT 
have significant observable consequences. 



The form of SCFT that has proved useful for the anal- 
ysis of experimental data is a phenomenological theory. 
It contains several parameters whose values are sensitive 
to details of monomer scale structure, which must be de- 
termined by comparison to experiment. In the simplest 
form of SCFT for systems with two types of monomers, 
these parameters are a Flory Huggins interaction param- 
eter x(T), and the statistical segment lengths of both 
species. SCFT does not predict how these parameters 
depend upon the details of molecular structure. Its use- 
fulness arises instead from its ability to predict phase be- 
havior, equilibrium structures, and diffuse scattering for 
systems containing polymers that are comprised of the 
same types of monomers, but that have different molec- 
ular weights and/or architectures. 

Any attempt to systematically calculate corrections 
to this SCFT, however, must start with some sort of 
micro-mechanical model. (We need a Hamiltonian to 
do statistical mechanics.) Here, as in several previous 
studies£&£&2ii& we start from a coarse-grained model 
of Gaussian chains with pairwise-additive interactions 
between monomers. Each coarse-grained monomer in 
such a model represents a subchain of many chemical 
repeat units, within a chain that contains many such 
monomers. Such models are thus implicitly coarse- 
grained to some cutoff length intermediate between the 
chemical monomer size and the polymer coil size. 

The long wavelength composition fluctuations that be- 
come important near the critical point of a blend, or the 
ODT of a symmetric diblock copolymer, exhibit a univer- 
sal phenomenology of their own. For a blend, sufficiently 
close to the critical point, this is the critical behavior of 
the Ising universality class. One might hope that our 
theoretical description of these long wavelength fluctua- 
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tions would be insensitive to the value chosen for a cutoff 
length, or to other arbitrary details of how our coarse- 
grained model behaves at very short length scales. Un- 
fortunately, this is not so, at least not in the simplest 
sense: Numerical values predicted for a variety of quan- 
tities turn out to be very sensitive to the value chosen for 
the cutoff length. The purpose of this paper is to show 
how this may be remedied by an appropriate renormal- 
ization scheme. 



A. Field Theory and Mean-Field Theory 

The calculations presented here make use of the auxil- 
iary field representation of the partition function that was 
introduced into polymer physics by Edward s 11 ^ 2 . This 
approach has been used in in several previous studies of 
fluctuation effects in polymer blends£&£>£. Let 

Z = J D\R]e~ u W (I) 

denote the partition function for a model with a poten- 
tial energy U, where J D\R] denotes an integral over all 
particle positions. The auxiliary field approach makes 
use of an exact transformation of the partition function 
for any model in which U is a sum of intramolecular po- 
tential and a pairwise additive potential for interactions 
between monomers. This transformation can be applied 
to either the canonical or grand-canonical partition func- 
tion. The transformation yields a representation of Z (in 
either ensemble) as a functional integral of the form 

Z = J D[J] e L [ J l (2) 

where J is an auxiliary field (or a pair of such fields, 
one for each monomer type) that has units of monomer 
chemical potential. This approach is discussed in more 
detail in Sec. [HO 

A saddle-point approximation to the auxiliary field 
functional integral is known to yield a very simple form of 
mean-field theory. The free energy functional obtained in 
this saddle-point approximation is the sum of the free en- 
ergy of an ideal gas of polymers plus a mean field approx- 
imation for the interaction energy. The average interac- 
tion energy obtained in this approximation is the same 
as that obtained by a "random mixing" approximation 
in which we ignore all correlations among the monomer 
positions. The underlying assumption of microscopically 
random mixing is analogous to that used in the Poisson- 
Boltzmann theory of electrolytes, or the original Flory- 
Huggins lattice model, which both thus "mean-field" the- 
ories in the same sense. Like other microscopic mean-field 
theories of this type, the resulting theory makes very sim- 
ple, but generally inaccurate, predictions about the rela- 
tionship between microscopic interaction parameters and 
macroscopic parameters, such as the effective x param- 
eter observed in scattering experiments. We will make 



a distinction in what follows between this microscopic 
mean-field theory and the phenomenological SCFT that 
is used to fit experimental data, which contains several 
adjustable parameters. 



B. UV Divergence and Renormalization 

One appealing feature of the auxiliary field approach 
is that it allows the effects of fluctuations of the aux- 
iliary field about this mean-field approximation to be 
treated by standard methods of diagrammatic pertur- 
bation theory, analogous to those used in the study of 
critical phenomenal Several studies have attempted to 
calculate corrections to the simple mean field theory for 
binary polymers blends^&i^ by introducing a Gaussian 
approximation for distribution of fluctuations of the aux- 
iliary J about its saddle point. In any perturbative field 
theory, a Gaussian approximation for fluctuations about 
the saddle point can be expressed diagrammatically in 
terms of Feynman diagrams that involve only a single 
"loop" , or a single wavevector integration. For this rea- 
son, this approximation is often referred to (and will be 
referred to here) as a "one-loop" approximation. 

In all of these calculations, it was found that the pre- 
dictions of the one-loop approximation for corrections to 
the mean-field free energy density and for the inverse 
structure factor^i^ are dominated by the contributions 
of short-wavelength fluctuations, with wavelengths of or- 
der the coarse-graining or (equivalently) a cutoff length. 
If the relevant Fourier integrals are cut off at a cutoff 
wavenumber A, predictions for the apparent x parame- 
ter are found to contain terms proportional to A and A 3 , 
which diverge in the limit A — > oo. In the jargon of field 
theory, the one-loop approximation was thus found to be 
ultraviolet (UV) divergent. 

This UV divergence is not peculiar to studies that rely 
on the auxiliary field approach, but arises in all coarse- 
grained models of fluctuation effects in polymer liquids. 
The most influential theory of fluctuation effects in di- 
block copolymer blends is that of Brazovskiii^., Fredrick- 
son and Helfand 4 (BFH). The BFH theory and its 
descendant s 15 ' 16 ! 17 ' 18 ? 19 are based on an effective Hamil- 
tonian formalism in which the partition function is ap- 
proximated as a functional integral with respect to a fluc- 
tuating monomer concentration field, rather than with 
respect to a fluctuating chemical potential. A saddle- 
point approximation for the functional integral used in 
the effective Hamiltonian approach yields the same mean- 
field theory as that obtained by a saddle-point approxi- 
mation in the auxiliary field approach. The BFH theory 
is based on a self-consistent one-loop approximations for 
fluctuations about this saddle-point. As emphasized by 
Kudlay and Stepanow-i£, this approach leads to UV di- 
vergences analogous to those encountered in the auxiliary 
field approach. 

The physical reason for this sensitivity to the value cho- 
sen for the cutoff length A -1 is not hard to understand. 
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The total free energy of a polymer liquid is only subtly 
different from that of a corresponding liquid of oligomers, 
because both are dominated by the effects of monomer 
scale liquid structure. It should thus be no surprise that 
the free energy density of a coarse-grained model is sen- 
sitive to our choice of cutoff length: The total free energy 
is sensitive to every detail of local fluid structure, at the 
shortest wavelength relevant to whatever model we con- 
sider. 

Phenomenological SCFT assumes that this sensitivity 
to local structure can be encapsulated within a few phe- 
nomenological parameters. Phenomenological SCFT is 
widely believed to be asymptotically exact in the high 
molecular- weight limit, except within a Ginzburg region 
very near the spinodal. If we asssume this to be true, 
however, we should nonetheless expect to obtain different 
estimates for the values of the SCFT phenomenological 
parameters from different mathematical approximations 
for the properties of a given microscopic model. Specif- 
ically, we should expect to obtain different approxima- 
tions for the x parameter and statistical segment lengths 
from a one-loop approximation than those obtained from 
the simple mean-field approximation. 

One-loop corrections to the predictions of the simplest 
mean-field theory should thus be understood to contain 
two conceptually different kinds of corrections: 

1. Corrections to mean- field estimates of the parame- 
ters of SCFT, e.g., of the statistical segment lengths 
and effective interaction parameters. 

2. Corrections to the phenomenology of SCFT, which 
cannot be absorbed into corrections to these pa- 
rameters. We expect these to become small in the 
limit N — > oo, except within a narrow temperature 
window near the spinodal. 

We expect corrections to the SCFT parameters to be sen- 
sitive to our treatment of short wavelength correlations, 
and thus, in a coarse-grained theory, to our choice of cut- 
off wavenumber. Conversely, we expect corrections to the 
phenomenology of SCFT, such as corrections to Gaussian 
chain statistics or to RPA predictions of the functional 
form of S{k), to arise primarily from longer wavelength 
fluctuations, and to be independent of A. One goal of 
this paper is to provide mathematical evidence for the 
consistency of this physical picture. 

Our ultimate goal is to construct a renormalized per- 
turbation theory that allows us to unambiguously calcu- 
late corrections to phenomenological SCFT. We assume 
(subject to confirmation) that a SCFT with renormalized 
parameters is asymptotically exact in the limit N — > oo. 
We thus hope to construct a theory in which all cor- 
rections to this form of SCFT can be shown to vanish 
in the limit N — > oo. The assumption that SCFT be- 
comes exact as N — > oo implies, however, that all large, 
UV divergent contributions to the calculated correlation 
functions (which generally do not vanish in the limit 
N — > oo) must preserve the dependence on wavevector, 



chain length, and architecture predicted by SCFT. This 
is possible only if all of the UV divergent terms in the 
one-loop approximation can somehow be absorbed into 
corrections to the values of the phenomenological param- 
eters used in SCFT. 

This criteria for "renormalizability" imposes some non- 
trivial constraints on the allowed functional forms of UV 
divergent corrections. It implies, for instance, that any 
UV divergent contribution to S^ 1 (k) that we wish to 
interpret as a renormalization of x be completely inde- 
pendent of k and N, and that it have the same value 
in a polymer blend and a diblock copolymer melt of 
the same composition, like the x parameter in SCFT. 
Similarly, it requires that any contribution to S (k) 
that we associate with a renormalization of a statisti- 
cal segment length exhibit the nontrivial but predictable 
/c-dependence of the change in the RPA prediction for 
S^ 1 (k) that would be caused by a slight swelling or a 
contraction of Gaussian chains due to a change in statis- 
tical segment length. We show here that these criteria 
are actually satisfied by the one-loop approximation. 



C. UV Divergences in Prior Work 

The UV divergence of the one loop theory has inspired 
a variety of responses in prior work. 

In studies of fluctuation effects in polymer blends, sev- 
eral authors have introduced a cutoff A that is assumed to 
be of order the inverse monomer size, and simply reported 
the dominant UV divergent contributions as functions of 
A.— i— . These results have sometimes been interpreted 
as meaningful predictions for the dependence of the x 
parameter^ii or the statistical segment length^ upon 
microscopic structure. We believe that this interpreta- 
tion is misguided: When a prediction of a coarse-grained 
model for some quantity depends upon a microscopic cut- 
off length, it indicates only that the value of that quantity 
is sensitive to details of local fluid structure that such a 
model should not be expected to accurately describe. 

Holyst and Vilgis^^ have instead argued for the in- 
troduction a cutoff length of order the polymer coil size 
R oc y/Nb as a way of controlling the UV divergence. 
Their reasoning is worth recounting: Holyst and Vilgis 9 
posed the question of whether the cutoff length that was 
needed as a result of the UV divergence of their theory 
should be taken to be a length of order the monomer 
size b or the coil size R. They noted that a cutoff length 
of order b would lead to corrections to the mean field 
(i.e., saddle-point) theory that do not become small in 
the limit N — > oo. They thus rejected this option in 
favor of the introduction of an ad hoc cutoff length of 
order R. These authors were thus led into a quandary, 
in part, because they did not make the conceptual dis- 
tinction made here between the microscopic mean field 
theory, which is certainly not exact in the limit N — > oo, 
and phenomenological SCFT theory, which (we argue) is. 

In Brazovskii's analysis of weakly first order crystal- 
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lization, he considers a UV divergent integral expression 
expression for S (k), but never mentions its divergence. 
This UV divergence is also not mentioned by Fredrickson 
and Helfand^ in their application of Brazovskii's analysis 
to diblock copolymer melts. Instead, these authors all 
report the UV-convergent part of the integral, which de- 
velops an infrared (IR) divergence at the SCFT spinodal, 
and discard the UV divergent contribution without com- 
ment. This approach is consistent with that normally 
taken in field theoretic studies of, e.g., the Ising critical 
point^i, in which it is well known that an analogous UV 
divergence of the one-loop theory can be absorbed into a 
renormalization of the critical temperature. This inter- 
pretation of the divergence in the Brazovskii model was 
made more explicit by Dobrynin and Erukhimovich 15 , 
who noted (in the appendix) it could also be absorbed 
into a shift in the critical temperature. 

This approach is the only appropriate one in the study 
of very generic models, such as the original Brazovskii 
model, in which the critical temperature is treated from 
the outset as an unknown parameter. It becomes prob- 
lematic only when such a field theoretic treatment of 
fluctuations is grafted onto a SCFT theory that we ex- 
pect to become exact in the limit N — > oo. SCFT pre- 
dicts nontrivial relationships between critical tempera- 
tures and order-disorder temperatures of "homologous" 
systems containing polymers that are constructed from 
the same monomers but that have different lengths or 
architectures (e.g., blends and diblock copolymer melts). 
If we were to treat the transition temperatures for dif- 
ferent such systems as completely independent parame- 
ters, we would forego the ability to say anything about 
corrections to SCFT predictions for these relationships. 
UV divergent contributions to the one-loop theory can 
be made truly benign only if they can be related to the 
parameters of SCFT. 

We are aware of two previous attempts to renormalize 
the theory of fluctuations in blends or copolymer melts, 
similar in spirit to that given here: 

Kudlay and Stepanow-i^ attempted to renormalize a 
refined version of the one-loop approximation for S^ 1 (k) 
in a diblock copolymer melt introduced by Fredrickson 
and Helfand. These authors proposed (in effect) that 
the UV divergence of S~ 1 (k) could be tamed if it could 
be absorbed into a renormalization of the Flory-Huggins 
X parameter. They concluded, however, that this inter- 
petation was not tenable for the theory that they consid- 
ered, because the UV divergence of S~ 1 (k) in this theory 
was shown to exhibit a nontrivial dependence on k, and 
because different results for the UV divergent contribu- 
tion were obtained for a diblock copolymer melt and for 
a binary blend of the same composition. 

The first successful attempt to renormalize a predic- 
tion of a one-loop theory was given by Wang. 8 Wang 
used the Edwards' auxiliary field approach to derive an 
expression for a one-loop correction to the k = limit of 
S~ 1 (k) in a binary homopolymer blend. He showed that 
the result was the sum of a UV divergent part that he 



interpreted as a renormalization of x-> and a UV conver- 
gent contribution that vanishes in the limit N — > oo, as 
suggested by the physical picture discussed above. The 
main limitation of Wang's calculation is that the method 
that he and others£&£ have used to calculate S^ 1 (k = 0) 
in a blend is not easily generalizable to allow the cal- 
culation of S^ 1 (k) at k 0, or to study fluctuations 
in block copolymer melts. Wang was thus not able to 
examine either the fc-dependent or the architecture de- 
pendence of his results. As such, Wang's calculation left 
open the question of whether his proposed renormaliza- 
tion scheme would have passed the more stringent con- 
sistency requirements imposed by Kudlay and Stepanow. 
Wang's analysis is discussed in detail in Sec. IIIIl Our 
analysis builds directly upon Wang's, and removes many 
of its limitations. 



D. Outline 

The paper is organized as follows: Sec. |TT] introduces 
the coarse-grained model of interest, as well as several 
mathematical and conceptual definitions that are needed 
to discuss our results. Sec. IIIIl contains an overview 
of our renormalization procedure and results, beginning 
with a review of Wang's results. We hope that readers 
who are interested in understanding only the main phys- 
ical ideas and results will be able to follow to this point. 

Sec. IIVIIXIII present the technical details of our anal- 
ysis. Sec. IIVI reviews the auxiliary field method, and 
the Gaussian/one-loop approximation. Sec. [V] reviews 
the one-loop calculation of the free energy density of a 
homogeneous binary blend, and a corresponding analysis 
of S' _1 (0), in which we retain some subdominant terms 
that previous authors ignored. In IVII we give a self- 
contained derivation of the one-loop correction to S (k) 
at arbitrary k by functional differentiation of the Gaus- 
sian approximation for the free energy functional. In 
IVIII we review the general diagrammatic rules obtained 
in Ref. [13]. These provide an alternate path to the 
same expressions for S(k), and also allow us to sepa- 
rate intramolecular from intermolecular correlations. In 
Sec. I Villi we convert our expression for the one-loop 
correction to S^ 1 (k) at fixed chemical potential (grand- 
canonical ensemble) to a corresponding correction for a 
closed system (canonical ensemble). In Sec. [Xj we an- 
alyze the UV divergence of the one-loop approximation 
for S~ 1 (k) for a binary blend. In Sec. IXI[ we present the 
corresponding calculation for a diblock copolymer melt. 
In Sec. IXTIl we show that 0(1/N) corrections to the UV 
divergent part of S~ 1 (k) can be consistently interpreted 
as the result of end effects, and of a square-gradient con- 
tribution to the interaction free energy. 

Sec. IXIIII presents a power counting analysis of the 
order of UV divergences of an arbitrary diagram, at any 
order in a loop expansion. We also show there that a 
renormalized loop expansion will yield an asymptotic ex- 
pansion of corrections to SCFT in powers of 1/y/N. Sec. 
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IXIVI discusses the relationship between this work and the 
extensive literature on excluded volume problem in poly- 
mer solutions. Concluding remarks are given in Sec. IXVI 



II. MODEL AND DEFINITIONS 

In what follows, we explicitly consider binary ho- 
mopolymer blends and diblock copolymer melts. We 
adopt a notation that allows for a mixture of any num- 
ber of molecular species constructed from a palette of any 
number of monomer types. Let p a be the number concen- 
tration of molecules of species a in a mixture. Let Ci(r) 
be the fluctuating number concentration of monomers of 
type i at point r. If H am i(s) denotes the position of 
monomer s of monomer type i on molecule number m of 
species a, then 



Ci(r) 



E 

sma 



<5(r - Rami(s)) 



(3) 



Here, the sum over s is taken over monomers of type i 
on molecules of species a, the sum over m is taken over 
molecules of species a, and the sum over a is taken over 
all species that contain i monomers. We will consider 
nearly incompressible liquid with an average volume v 
per monomer, in which v is the same for all monomer 
types, and independent of composition in a mixture. 



A. A Coarse-Grained Model 

Consider a coarse-grained model for polymer liquids in 
which the total potential energy is a sum 



U = [/chain + CAnt + [Text [h] , 



(4) 



of an intramolecular potential energy f/ c ham, a pair inter- 
action potential 

Uint = \ Y,j dr j dr ' U ^ V - r 'Mr) Ci (r') (5) 



and an external potential 



E 



dr hi(r)ci(r) 



(6) 



Here, J7y(r — r') is a pair potential for interactions be- 
tween monomers of types i and j, and hi(r) is an external 
potential field conjugate to Cj(r). The external potential 
fields are introduced as a mathematical convenience, and 
are set to zero in all final expressions. 

In what follows, we sometimes use a compact nota- 
tion in which a binary operator '*' is used to indicate 
integration over a shared coordinate (or wavevector) and 
summation over a shared monomer type index. For ex- 
ample, 



Ui n t = -c*U *c 



U 



ext 



-h * C 



(7) 
(8) 



in this notation. 

We assume in what follows that C/ C hain is adequately 
approximated at the length scales of interest by the 
stretching energy of a continuous Gaussian chain. Let 
hi denote the statistical segment for monomers of type i 
on an isolated chain (with ?7i n t = 0), and 



k = v/bf 



(9) 



be the corresponding packing length for monomers of 
type i. 

We consider a class of models in which the Fourier 
transform Uij (is.) of the pair potential ;(r — r') is of the 
form 



U ij (k) = U ij F(k/A) 



(10) 



where Uij is a matrix of interaction strengths with di- 
mensions of volume, A is an inverse range of interaction, 
and F(k/A) is a function that approaches 1 as k/A — > 0. 
This corresponds to a pair potential in coordinate space 



£Z y -(r-r') = £7yA 3 F(A|r - r' 



(11) 



where the function F is the inverse Fourier transform of 
F, and satisfies normalization condition J c?x F(x) = 1. 
The interaction matrix Ua is taken to be of the form 



Uij = v 



Bo Bo + Xo 
Bo + Xo Bo 



(12) 



The resulting potential energy reduces in the limit of slow 
spatial variations, in which the characteristic scale for 
gradients is much larger than A -1 , to a continuum ap- 
proximation 



Un 



dr 1 t: S o(ci + c 2 ) 



X0C1C2 



(13) 



in which Bo/v is the mean-field compression modulus. 

The analysis presented here can be carried out in either 
canonical or grand-canonical ensemble, with only minor 
differences. We will work in grand-canonical ensemble. 
Let Z[h] denote the grand canonical partition function, 
for some choice of chemical potentials. Let -F[(c)] be the 
corresponding free energy expressed as a functional of the 
average monomer concentration fields. This is defined by 
the Legendre transform 



F[(c)] =- In Z[h] + h*(c) 



(14) 



Here, h is the external field required to establish a 
monomer concentration (c), which is related to F by a 
functional derivative 



hi(r) 



SF[(c)} 



(15) 



In Eq. (|14p . and hereafter, we use energy units in which 
k B T = 1. 
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B. Correlation Functions 

We focus on the calculation of the correlation function 
S y (r,r') = (^(r)^(r')) (16) 
and its Fourier transform 



Sy(k) = J dv 1 S l3 {r,0y* r , 



(17) 



where Sa(r) = Cj(r) — (cj(r)). This quantity obeys the 
identities 



8 2 hxZ[h] 
Shi^Shjir') 
S 2 F[(c)} 
5(c i (r))6{c j (T')) 



(18) 
(19) 



The inverse S^ 1 is defined in coordinate space by requir- 
ing that S * S 1-1 = 8, where 8 denotes 6(r,r')8ik, or in 
the Fourier space by requiring (for a homogenous liquid) 
thatE i 5y 1 (k)5 ifc (k)=^ fc . 

We define an intramolecular correlation 



Q a jj(r,r') = ^(c ami (r) c amj (r')) , 



(20) 



that describes correlations between pairs of monomers on 
the same molecule to of a specified species a, in which 



(r)=^£(r-R ami (s)) 



(21) 



is the concentration of monomers of type i on a specific 
molecule to of species a. The sum over to in Eq. ([20]) is 
over all molecules of type a. The sum over s in Eq. Q21j) 
is over all monomers of type i on molecule to. 

In a molecular mixture, let f2jj(r, r'), with no species 
index a. denote the total contribution 



^•(r,r') 



E 



n a ,y(r,r') , 



(22) 



of intramolecular correlations to Sy(r, r'). Here, the 
sum is taken over all species a that contain both i and 
j monomers. In a binary mixture of homopolymers of 
species a = 1 and 2, the only nonzero element of for 
homopolymer 1 is the element i = j = 1, while the only 
nonzero elements of J7y are the diagonal elements, with 
i = j. In a single component copolymer melt, £l a ,ij and 
Ojj are identical. 



C. Self-Consistent Field Theory 

Polymer SCFT is a density functional theory of inho- 
mogeneous polymer liquids. It is based upon an approx- 
imation of the free energy functional -F[(c)] as a sum 



F[(c)} =F chain [(c)]+F int [(c)] 



(23) 



in which -Fchain is the free energy of a hypothetical refer- 
ence system of non-interacting polymers with a specified 
average concentration profile, and -Pint [(c)] is an addi- 
tional "interaction" free energy functional. 

The functional -F c hain[(c)] is the free energy of a hypo- 
thetical system of non-interacting chains, with U- ln t = 0, 
in which a monomer concentration (c) is maintained by 
a potential U cx t[h], with an applied field h. The field h 
thus satisfies an identity 



h(r) = 



8F 



chain 



[(c)] 



(24) 



as a special case of Eq. (|T5|) . Applying Eq. (fT5|) to Eq. 
([23]) for F[(c}} yields a self-consistent field equation 



hi(r) = hi(r) 



8F int [(c)} 



(25) 



Here, hi(r) is the external field that must be applied to 
the interacting fluid to establish a monomer concentra- 
tion field (c), hi(r) is the field required to establish the 
same concentration field in the non-interacting reference 
system, and 8Fi nt /8(c} is an "internal" field contribution. 

By itself, Eq. ([23| is merely a definition of i*l n t, 
rather than a predictive theory. The standard coarse- 
grained model for dense multi-component liquids of flex- 
ible polymers^ may be obtained by supplementing this 
with the following physical assumptions: 

a) Gaussian Chains: Polymer conformations are ade- 
quately described at the mesoscopic scales of inter- 
est as Gaussian random walks. 

b) Locality: The interaction free energy Pint may be 
adequately approximated by a local functional, of 
the form 



= J dr /i„t(ci(r),c 2 (r)) , 



(26) 



where /mt(ci(r), c 2 (r)) is a free energy density at 
point r that depends only upon the monomer con- 
centrations at point r. More precisely, it is assumed 
that the range of any nonlocality in _Fi n t is of order 
the monomer size, and that this is small compared 
to the length scales of interest in applications of the 
coarse-grained theory. 

c) Incompressibility: It is often assumed that liquid 
is essentially incompressible at the length scales of 
interest. If coarse-grained 1 and 2 monomers have 
been defined so as to occupy the same volume v, 
this may be summarized by a constraint 



= < Cl (r)) + (c 2 (r)) 



(27) 



The simplest form of phenomenological SCFT for sys- 
tems with two types of monomers assumes, in ad- 
dition, that /; nt in an incompressible liquid may be 
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adequately described by the Flory-Huggins expression, 
/int = v x( c A( r )}(cB( r )) with an interaction parameter 
X- This simple assumed composition dependence is not 
a necessary or essential part of the theory. 

To precisely define the decomposition of F into -Fchain 
and -Pint in Eq. (|23|) . we must make a choice of what 
single-chain reference Hamiltonian to use to define the 
ideal gas free energy Fchain- The physical reasoning un- 
derlying SCFT suggests that this reference Hamiltonian 
should yield chain conformations that, in a homogenous 
state, are similar to those in the liquid of interest. Even 
chains that are approximately Gaussian in a dense liq- 
uid generally have statistical segment lengths that are 
slightly different from those of a corresponding system 
of non-interacting chains.— In what follows, we will thus 
define F cna m to be the free energy of a hypothetical sys- 
tem of non-interacting chains in which the single-chain 
reference Hamiltonian is chosen so as to yield exactly the 
same single-chain probability distribution as that found 
for chains in the homogeneous state of the liquid of in- 
terest. The intramolecular reference Hamiltonian used 
to calculate F cna in should thus be understood to be a 
single-chain potential of mean force, rather than the bare 
intramolecular potential J7 c hain- Our results indicate that 
the UV divergent part of the one-loop contribution to the 
remaining free energy F- lrA is actually a local functional 
only if F c h a in and Fj nt are defined in this way. 



D. Ornstein-Zernicke Relation 

It is useful to introduce a generalized Ornstein- 
Zernicke (OZ) relation^^^ 

^ 1 (r,r')=Or>( r ,r')-a i (r ) r') (28) 



or 



^ 1 (k)=fi-. 1 (k)-C y (k) (29) 



in a homogeneous liquid. Eq. (|2"8")) defines the direct 
correlation function Cy . 

The single chain correlation function fly (k) is related 
to the reference free energy F c hain[(c)], as defined above, 
by an identity 



S 2 F chain [(c)} 
6( Ci (r))5{ Cj (r>)) 



(30) 



analogous to Eq (|19p . Differentiation of the free energy of 
a reference system of non- interacting molecules yields the 
intramolecular correlation function because the only cor- 
relations in such an ideal gas are intramolecular. Differ- 
entiation yields the true intramolecular correlation func- 
tion fir. (r, r') in the dense liquid of interest, rather than 
that of a gas of molecules with the bare intramolecular 
potential t/ c hain, as a result of the definition of F c h a in 
discussed above. 

The direct correlation function is intimately related 
to the interaction free energy Fiat defined by Eq. (f2"3")l . 



Combining Eqs. (JT9J and ([30) with Eq. (J23J) for F yields 
an identity 



<%(r,r') = - 



S 2 F int [{c)} 
S( Cl (r))d( Cl (r')) 



(31) 



Assuming that F nrt is a local functional thus implies that 
the Fourier transform dj (k) should be independent of 
k. 

In the case of a binary blend of two homopolymers, the 
k — > limit of the transform dj (k) is related to the the 
composition dependence of the free energy density of a 
homogeneous mixture. Consider a homogeneous blend of 
two homopolymers, in which Ni is the degree of polymer- 
ization of i homopolymers, Cj = NiPi is the macroscopic 
(i.e., spatial average) concentration of i monomers, and 
/int(ci,C2) is the interaction free energy density of the 
mixture. Then 



limC^k) 



dcidc-j 



(32) 



E. Incompressible Limit 

A simplified expression for Sij (k) may be obtained in 
the limit of an incompressible liquid. In a nearly in- 
compressible liquid containing two types of monomer of 
equal volume, the 2x2 matrix <SVj(k) has two widely 
disparate eigenvalues: In this limit, one eigenvector of 
SVy(k) must approach a pure "composition" fluctuation 
mode, [<5(ci(k)), <5(c 2 (k))] cx [1,-1], which satisfies the 
constraint <$(ci(k)) + <5(c 2 (k)) = 0, and thus avoids the 
large free energy penalty for changes in total monomer 
density. The other eigenvector must have a vanishing in- 
ner product with the first, and so must approach a pure 
"compression" mode, [<5(ci(k)), <5(c 2 (k))] cx [1,1]. 

The incompressible limit of Sij (k) may be obtained by 
assuming that the eigenvectors approach the limits de- 
scribed above, and taking the eigenvalue of the compres- 
sion mode to vanish. Alternatively, it may be obtained by 
assuming dj (k) to be of the form dj (k) = —B+SCij (k) , 
and taking B to infinity while keeping SCij (k) finite. Ei- 
ther method yields a matrix correlation function of the 
form 



S(k) = 5(k) 



+ 1 
-1 



with a scalar correlation function 

|0(k)| 



5(k) 



r> + (k)-2 UXa (k)|fi(k)| 



(33) 



(34) 



where 



n+(k) = fin(k) + fi 22 (k) + fi 12 (k) + r! 21 (k) 
|Q(k)| = nii(k)n 2a (k)-n 12 (k)n a i(k) (35) 
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are the sum of elements of fly(k) and its determinant, 
respectively, and where 

Xa (k) = ^[C u (k) + C 22 (k) - 2C 12 (k)] (36) 

is a wavenumber-dependent "apparent" x parame- 
ter. This definition was introduced by Schweizer and 
Curro 25 i 26 , using similar reasoning. 

In an incompressible binary homopolymer blend, the 
long- wavelength limit 



X a (0) = lim Xa(k) 

may be expressed as a derivative 

v <9 2 / int (>i) 



Xa(0) = -- 



(37) 



(38) 



where /i nt is expressed as a function of the total volume 
fraction <pi = vci of one of the two homopolymers. 



III. OVERVIEW 

In this section, we provide an overview of our renor- 
malization procedure, and summarize our main results. 



A. Mean Field Theory 

A saddle-point approximation to the Edwards auxil- 
iary field theory yields a very simple mean field the- 
ory. This is a form of SCFT in which -F c hain[(c)] is ap- 
proximated by the free energy functional for a gas of 
non-interacting chains, with a single-chain Hamiltonian 
C4hain, and in which 



F h 



dr dv 1 ( Ci (r))[/y(r-r')( Cj (r')) 



(39) 



This approximation yields a self-consistent (saddle-point) 
field 



hi(r) = hi{r) - / dr' U l3 {v - r')( Cj (r')) 



(40) 



In this approximation, the internal contribution to hi (r) 
(i.e., the convolution integral in the above) is liter- 
ally a "mean" field, insofar as it is approximated by 
the ensemble average of the fluctuating potential field 
— / dr'Uij(r — r')cj(r') at point r. 

The characteristic features of this microscopic mean- 
field theory (however it is obtained) are that: i) Fi nt is ap- 
proximated by an expression for the average interaction 
energy (£/i n t) that neglects all correlations in monomer 
density, and ii) intramolecular correlations are taken to 
be identical to those of a gas of non-interacting molecules. 
This theory makes very simple predictions about the 
relationship between SCFT parameters and the micro- 
scopic parameters: It yields a direct correlation function 
Cij(k) = —Uij(k), an apparent x-parameter Xa{k) — xo 
for k <C A, and statistical segment lengths equal to those 
of the non-interacting chains. 



B. One- loop Approximation for Free Energy 

The treatment of fluctuation effects given in this pa- 
per is based on a one-loop (or Gaussian) approximation. 
Several author s££££ have previously obtained a one-loop 
approximation for Xa(0) by calculating the free energy of 
a homogeneous blend as a function of composition, and 
applying Eq. (f3"5| . All of these studies started from an 
idealized model in which chains are treated as continuous 
Gaussian threads, the pair interaction is taken to be nom- 
inally point-like, as in the Edwards model of excluded 
volume interaction in solution, and in which the liquid 
is taken to be nominally incompressible. When these as- 
sumptions are all taken literally, the resulting one-loop 
correction to Xa(0) is given by a Fourier integral that di- 
verges at large wavenumber, i.e., that is UV divergent. 
The divergence can be removed either by introducing a 
discrete chain model or by introducing a nonzero range 
for two-body interactions*^ In all of these previous stud- 
ies, the integral expression for Xa(0) was regularized by 
restricting the integral to wavenumbers less than a cutoff 
wavenumber A. 

This regularized one-loop approximation for binary ho- 
mopolymer blends, which is presented in Sec. [Vj yields 
an apparent x parameter of the form 

2 TT 

Xa(0) = Xo + AA 3 + B X oA + J2w A + 6x *^ (41) 



where 



.4 
B 
Hi 
H 2 



(h-h) 2 v 
24n 2 l 2 

7T 2 P 



3hh_ 
2tt 2 P 

2n 2 T 3 



and 



l(4>i) = fail + fah 



(42) 



(43) 



Here, <5x*(0) is a contribution that remains finite in the 
limit A — > oo (i.e., that is UV convergent). 

The dominant UV divergent parts of this result have 
been reported previously. De la Cruz et al~ consid- 
ered the structurally symmetric case &i = £>2 7 for which 
A = 0, and reported the contribution of the form BxA. 
Fredrickson and Liu&£ instead considered the athermal 
case bi ^ b 2 and \o — 0, and so found the AA 3 con- 
tribution. Wang& retained both of these contributions. 
The terms proportional to A/A have not been retained 
or analyzed in previous work. 

Wang was the only one to attempt to renormalizc this 
theory, in the sense proposed here, by absorbing the 
strongly cutoff dependent contributions to Xa into a re- 
definition of the SCFT x parameter. Wang proposed that 
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the x parameter used in phenomenological SCFT theory 
be identified, within the one-loop approximation, with a 
sum 

Xcff = Xo + AA 3 + B X A . (44) 

This definition of XcS is independent of chain length and 
architecture, and is thus consistent with the physical pic- 
ture of x m SCFT as a parameter that is sensitive to 
details of local fluid structure, but insensitive to changes 
in chain connectivity at longer length scales. 

By absorbing the dominant cutoff-dependent terms 
into a renormalization of \, while ignoring the contribu- 
tions of order A/N, Wang was able to isolate the remain- 
ing UV convergent contribution #x*(0). It is this quan- 
tity that contains predictions of the renormalized theory 
for corrections to SCFT phenomenology. By introducing 
several further approximations, Wang obtained an ap- 
proximate analytic expression for <5x*(0). Both Wang's 
analytic approximation <5x*(0) and the full one-loop in- 
tegral expression for this quantity may be expressed in 
the non-dimcnsionalizcd form 

N6 X *(0) = j^5x*(xN,<t>uNi/N 2 ,b 1 /b 2 ) , (45) 

in which Sx* is a dimensionless function of all of the 
dimensionless variables relevant to SCFT. Here, N is a 
reference degree of polymerization (e.g., typically Ni or 
N 2 ), 

N = Nb 6 /v 2 (46) 

is an invariant degree of polymerization, and b is a ref- 
erence statistical segment length (e.g., b\ or 6 2 ). The 
value of x used in the integral that defines the RHS of 
Eq. (|45p may be taken to be either XcS, to obtain a 
renormalized perturbation theory, or may be replaced by 
Xa(0) to obtain the type of self-consistent one- loop ap- 
proximation considered by Wang. In either variant of the 
theory, the resulting correction to SCFT vanishes in the 
limit N — * oo as a result of the prefactor of A -1 / 2 . 

After showing how to remove the UV divergence from 
this calculation, Wang focused primarily on a discussion 
of 5x* (0) near the spinodal, which determines the critical 
behavior of a blend. We will defer our own discussion of 
this subject, and other physical predictions of the one- 
loop theory, to a subsequent publication. 

C. One-loop Approximation for S(k) 

In Sees. IXland lXll we analyze the UV divergent con- 
tributions to S (k) in homopolymer blends and diblock 
copolymers, respectively. The physical assumptions un- 
derlying phenomenological SCFT allow for the possibil- 
ity that both the x parameter and the statistical lengths 
may be different in a one-loop approximation from that in 
mean-field theory. Renormalization of a local x param- 
eter is expected to introduce a fc-independent change in 



S~ 1 (k). Changes in statistical segment length would give 
rise to changes in S~ 1 (k) with a non-trivial, but foresee- 
able, wavenumber dependence. To distinguish these ef- 
fects, we use the OZ expression for (k), and calculate 
separate one-loop contributions to the single chain corre- 
lation function f2y (fc) and to the direct correlation func- 
tion Cij{k). As in the simpler calculation of Xa(k = 0) 
discussed above, we divide the one-loop contribution to 
each of these quantities into a UV-divergent part that we 
calculate explicitly, and a remaining UV-convergent part 
that will be examined elsewhere. 



1. Intramolecular Correlations 

The standard SCFT assumes that polymers are ap- 
proximately Gaussian in a dense mixture, but not that 
the statistical segment lengths are necessarily the same 
as those in a reference system of non- interacting chains. 
Let flij (k; b) denote the single chain correlation function 
for a Gaussian chain with a specified statistical segment 
length b, or (for a diblock copolymer) with a specified pair 
of statistical segment, b = (61,62)- For a homopolymer, 
ttij(k) with i = j is proportional to a Debye function. If 
the only effect of interactions upon the single chain cor- 
relation function Qij(k) were to change the values of the 
statistical segment lengths, we would expect to find 

Q ij (fc)~Q ij (A;;6o + 56) , (47) 

where 60 is the "bare" statistical segment length (or 
lengths) for non-interacting chains, and Sb is a correction 
arising from interactions. This correction is calculated 
here to first order in a loop expansion. To first order 
in an expansion in powers of (56, or to first order in a 
loop expansion, we would expect a renormalization of 6 
to yield a result of the form 

Qijik) ~ %(fc; 6 ) + ; ^ Sb >* ■ ( 48 ) 

k 

In the case of a homopolymer blend, the only nonzero 
terms in the sum are those with i = j = k. The more 
general notation is required for a diblock copolymer melt, 
in which O12 (fc; 6) is a nonzero function of both 61 and 
62. In addition to this renormalization of the statistical 
length, we expect to find small corrections to Gaussian 
chain statistics. We thus expect to find a total one-loop 
contribution to 0^- (k) of the form 

sn ij{k) ~ J2 ™f£ b) 6b k + sn* 3 (k) (49) 

k 

in which dOjy (k) is a UV convergent correction that van- 
ishes in the limit N — ► 00. 

Our results for the one-loop contribution to fiy (k) are 
completely consistent with the above discussion: We find 
that the UV-divergent parts of the one-loop contribu- 
tion to flij(k) in both homopolymer blends and diblock 
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copolymer melts have precisely the wavenumber depen- 
dence suggested by Eq. (|4*5|) . The calculated fractional 
change in statistical segment length 6j, to first order in a 
loop expansion, is given by an expression 



5h 
bi 



n 2 l 



A 



(50) 



that increases linearly with A. Identical expressions for 
Sbk are obtained from calculations of Oy (fc) in a ho- 
mopolymer blend and in a disordered diblock copolymer 
melt with the same overall composition. Eq. ([517)) for 
Sbk was obtained previously by Wang2£ by considering 
the effect of fluctuations upon the end-to-end vector of a 
single chain in a binary blend. 

The UV convergent one-loop contribution SVt*j can be 
isolated by subtracting the above UV-divergence from 
the total one-loop contribution. Our result for this quan- 
tity is given by a convergent Fourier integral that can be 
expressed, in either blends or diblock copolymer melts, 
in non-dimensional form 



(k) = j^6ti* 3 (kR, X N,...,) (51) 



in which SO^a is a dimensionless function of all of the di- 
mensionless variables of SCFT (i.e., x_/V, <f>, Ni/N 2 , and 
61/62 in a blend, or of x^ \ J ■> an d 61/62 in a diblock 
copolymer melt), and of a non-dimensionalized wavevec- 
tor kR. Here, R = V~Nb, N is a reference degree of poly- 
merization (i.e., Ni or N2 in a blend or TV in a diblock 
copolymer melt), and 6 is a reference statistical segment 
length (i.e., 61 or 6 2 ). 



2. Inter-molecular Correlations 

We now consider inter-molecular correlations, as char- 
acterized by the direct correlation function Cy(fc). We 
showed in subsection III Dl that Cy (k) can be expressed 
as a functional derivative of the SCFT interaction free 
energy F mt . The standard SCFT assumes that Fi nt is a 
local functional of monomer concentration fields, and is 
independent of chain length and architecture. If these as- 
sumptions were rigorously correct, the quantities Cij (k) 
and Xa(k) would be completely independent of k, chain 
length and architecture, but could depend upon compo- 
sition. 

In Sec. [Xj we show that the one-loop approximation 
for Xa(k) in a binary blend can be expressed as a sum 



Xa(k) =XcS + Y^ 



HjjkRi) 



A + S x *(k) (52) 



in which XeS is the renormalized x parameter given ex- 
plicitly in Eq. ([44]) . which is independent of k. Here, 
Sx*{k) is a UV convergent contribution that becomes sig- 
nificant only near a spinodal, while Hi(kRi) is a dimen- 
sionless function that approaches the constant Hi given 



in Eq. (|4ip in the limit k — > 0. The corresponding ex- 
pression for Xa{k) in a diblock copolymer melt, which is 
analyzed in Sec. IXI1 is of the form 



Xa{k) = Xcff + 



H(kR) 
N 



A + S X *(k) 



(53) 



where TV is the length of the diblock. Both calculations 
confirm that the contributions to S^ 1 (k = 0) that Wang 
absorbed into Xeff do indeed correspond, for k ^ 0, to 
quantities that are independent of k and chain architec- 
ture, consistent with his physical interpretation. 

Our results for the UV divergent part of Xa{k) simplify 
considerably in the limit £1 = 1% = I of equal statistical 
segment lengths. In this case, the coefficient A vanishes 
in Eq. (144)) for XcS- In addition, we find that the co- 
efficients Hi (in a binary blend) and H (in a diblock 
copolymer melt) also vanish in this limit, for all k. In 
this case, we thus obtain 



Xa(k) = x 



1 - ° IA 

7T J 



<*X*(k) , 



(54) 



where I = l\ = 1%. Here, <5x*(k) is a UV convergent 
part. The term proportional to A is the UV divergent 
correction to x that was originally identified by de la 
Cruz et al^. 

The corrections to SCFT that are of physical inter- 
est arise from the smaller UV-convergent contribution 
Sx*(k). The exact one- loop expression for this quantity 
is given by a UV convergent integral, that may be ex- 
pressed, in either blends or diblock copolymer melts, in 
the non-dimensional form 



N6 X * a (k) 



1 



A71/2 



8x*{kR, x N,...,) 



(55) 



where Sx* is a dimensionless function of the same vari- 
ables as those used in Eq. (fSTj) . The k — * limit of our 
result for this quantity agrees with that obtained by ex- 
amining the composition dependence of the free energy 
density. 



D. End-Effects and Square-Gradient Terms 

The UV divergent terms in Eqs. (|52|) and (|53|) that are 
of order A/N cannot be absorbed into the renormalized 
SCFT parameter XeS, because they depend on molecular 
weight N, wavenumber k, and chain architecture. To 
explain the physical origin of these terms, we consider a 
slightly generalized form of SCFT in which we allow for 
two physical effects that are not included in the standard 
form of the theory. These are: 

1. Excess free energies associated with chain ends and 
junctions between blocks in block copoymers. 

2. A square-gradient contribution to the one-loop in- 
teraction free energy. 
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In Sec. IXII1 we thus compare our one loop results for 
Xa{k) to the predictions of a phenomcnological model in 
which the UV divergent part of the one-loop contribution 
to -F^t is assumed to be of the form 



SF int = / dr[ Sf 



h.,*\ + ^ J d a ilj a + -D{V(t>if) , (56) 



where 4>i(r) = vci(v) is a local volume fraction for one 
of the components in an incompressible liquid. Here, 
^/locai(r) is a local free energy density, d a (r) is a lo- 
cal concentration of chain end or junction 'defects', and 
tpa(r) is an excess free energy arising from the presence 
of a 'defect' of type a. In a binary homopolymer blend 
in a binary homopolymer blend, d\ (r) and di (r) are the 
concentrations of chain ends for chains of type 1 and 2, 
respectively. In a diblock copolymer melt, the index a 
can take values 1, 2, or J, where dj(r) is the local concen- 
tration of junctions between the blocks. The free energy 
density £/i oca i(r), the defect free energy ?/> a (r) for each 
type of defect, and the coefficient D are all assumed to 
be sensitive to local fluid structure, and so may depend 
upon the cutoff A, the statistical segment lengths, and 
the local composition </>i(r). 

Excess free energies for chain ends and junctions can 
arise in the one-loop approximation, even in a model in 
which the end and junction monomers are assumed to 
be identical to other monomers of the same type, simply 
because the local environment of a chain end or junction 
is different from that of a monomer in the middle of a 
long chain. This difference is captured at a crude level 
even by a one-loop theory. In any real polymer liquid, 
the excess free energy associated with chain ends would 
also be sensitive to any differences between the actual 
chemical structure of the terminal units and the chemical 
repeat unit. 

Previous analyses of the one-loop approximation have 
ignored these 0(A/N) contributions to x a (k). It was 
tempting for us to do the same, on the grounds that these 
terms are smaller by a factor of 1 /N than those absorbed 
into Xcff- Actually, however, these contributions are the 
same order in an expansion in l/N as the term propor- 
tional to xo-, since xo must be less than a binodal value of 
of (D(l/N) in order for the homogeneous state of interest 
to remain stable. We also found that we needed to ana- 
lyze and subtract these 0{A/N) divergences in order to 
carry out the numerical renormalization procedure that 
we are now using to calculate UV- convergent predictions 
of the theory. This procedure is described briefly in Sec. 

EH 

We discuss contributions to Sx a (k) that arise from a 
postulated square-gradient contribution to F- mt simulta- 
neously with those that arise from end and junction de- 
fects because we find that these contributions are other- 
wise difficult to disentangle. The square-gradient contri- 
bution to 5Fi nt in the above model simply adds a contri- 
bution Sxa(k) = — %Dk 2 to Xa{k). This, however, can 
also be written as -§£> (kR) 2 /Nb 2 , where R = \/Nb, 



and so can be absorbed into a contribution of the more 
general functional form H(kR)A/N, if D cx A. The only 
way for us to identify a squared-gradient contribution 
is thus to explicitly calculate the wavenumber depen- 
dence that would be produced by end and junction de- 
fects alone, and then see if our one loop results for Sxa(k) 
can be expressed as the sum of this defects contribution 
plus an additional square gradient contribution. 



1. Free Energy of Homogeneous Liquids 

Expressions for <5/i OC ai and i/j a in this generalized SCFT 
may be obtained by examining the one-loop contribution 
to the free energy density of a homogeneous liquid. 

The quantity £/i OC ai may be obtained by considering 
the N — > co limit of the UV divergent part of the free 
energy density of either a homogeneous binary blends or 
a disordered diblock copolymer melts. In either case, we 
obtain 



Sf\o 



1 



cal 



12tt 3 
ir 2 v I 



In 12,Bo 



I 

vA 2 



A a 



—n —XO^l&A 



(57) 



where <pi is the macroscopic volume fraction of i 
monomers. The one-loop contribution to the parameter 
Xeff defined in Eq. |44|) is related to <5/i OC ai by a second 
derivative with respect to <j>\, as in Eq. (|38[) . 

In a homogeneous binary blend, we find that the total 
UV divergent part of the one-loop free energy density Sf 
of a liquid of finite chains may be expressed as a sum 



Sf = Sf local + 



20i 
Niv 



202 

N 2 v 



4>2 



(58) 



where 2<f> a /N a v is the concentration of chain ends for 
chains of type a, and 



i> 



4tt 2 1 



llA 



(59) 



The quantity ip a is thus tentatively identified as the ex- 
cess free energy of an a chain end. 

The corresponding free energy for a disordered diblock 
copolymer melt can be written as a sum 

Sf = 5/local + -J- (ipl + ^2 + i/>j) (60) 

Nv 

where ip\ and ^2 are excess free energies for the ends of 
the 1 and 2 blocks, as given by Eq. (|BT)|) . The remaining 
energy ipj is given by 



At: 2 I 



(h - kfA 



(61) 



This is assumed to be the excess free energy arising from 
the junction in a single diblock copolymer. The fact that 
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this expression for ipj vanishes in the limit l± = I2 is con- 
sistent with the fact that in the limit Xo = an d h = h 
the two blocks of a diblock copolymer become indistin- 
guishable, so that the excess free energy associated with 
the junction must thus vanish in this limit. (This argu- 
ment would allow ipj 7^ for l\ — I2 and xo 7^ 0, but this 
does occur to first order in a loop expansion.) 

2. Composition Fluctuations 



Because both the defect contributions and this square- 
gradient coefficient vanish when li = l 2 , the 0(A/N) 
contribution to JS 1-1 ^) vanishes in this case. 

The success of this approach strongly suggests that the 
UV divergent part of the one-loop contribution to the free 
energy functional -Pint [(c)] is of the form assumed in Eq. 
(I56p . even for strongly inhomogeneous liquids. We thus 
hope that this result will also provide a basis for removing 
the UV divergence of the one-loop approximation for the 
free energy of ordered phases of block copolymer melts. 



The main evidence for our interpretation of the 
0(A/N) contributions to S~ 1 (k) as a combination of end 
and junction effects and square-gradient contributions is 
a demonstration that the wavenumber and parameter de- 
pendence of these terms can be explained by this inter- 
pretation. 

In Sec lXIII we present an RPA calculation for S~ 1 (k) in 
a system with an additional free energy of the form given 
in Eq. ([56]) . using the explicit expressions for ip a given in 
Eqs. (|59|) and (f6Tj) . For a binary blend, the required cal- 
culation is the same as the RPA calculation for a blend of 
polymers with chemically distinct end-groups, or a mix- 
ture of CAC and DBD triblocks with very short C and 
D end groups. The calculation for diblock copolymers 
is similar to that for a pentablock copolymer with very 
short end and middle blocks. The inclusion of end effects 
in the RPA generally yields contributions to S~ 1 (k) that 
are proportional to 1/N, with a nontrivial ^-dependence 
that is different in binary blends and diblock copolymer 
melts. 

The excess free energies associated with chain ends and 
junctions affect the collective correlation function S(k) if 
and only if the defect free energies ip a (r) depends upon 
the composition 4>i(r) of its immediate environment. If 
the free energy per defect depends upon the composition 
of its environment, defects will tend to cluster where their 
free energies are lowest, and to favor collective composi- 
tion fluctuations that lower the total defect free energy. 
In systems with l\ ^ I2, the above expressions for ip a (r) 
generally do depend upon </>i(r), because of the compo- 
sition dependence of the average packing length l(4>i). 
Both end and junction defect free energies are indepen- 
dent of 4>\ when l\ = I2, implying (correctly) that we 
should find no "defect" corrections to S~ 1 (k) in this spe- 
cial case. 

We find that the results of this generalized RPA cal- 
culation exactly reproduce the structure of the UV- 
divergent 0(1/ N) terms obtained in our one-loop cal- 
culation of S^ 1 (k), if we allow for both end effects and a 
square-gradient contribution to .Pint [(c)]. The contribu- 
tions to S~ 1 (k) that arise from the defects vanish in the 
limit l\ = I2, for the reason discussed above. The value 
of the coefficient D that we infer by this method (which 
is the same in binary blend and diblock copolymer melts 
of the same composition) is given by 



D = 



{h-hf 
3tt 2 / 2 



A 



(62) 



IV. AUXILIARY FIELD METHOD 

The Edwards' functional integral representation of Z 
may be obtained from the identity 



e -h c * u * c = TV 



,-i J*U~ 1 *J+iJ*c 



(63) 



Here, J D[J] represents a functional integral with respect 
to an auxiliary chemical potential field, where Jj(r) is a 
field component that couples to Cj(r). The constant N is 
given by the integral 

A = f D[J] e-h^U-^J ( 64 ) 



Substituting this representation of e~ Uint into the defini- 
tion of the grand-canonical partition function Z yields a 
functional integral 



Z[h] = N- 1 Jd[J] e L[h - J] 



(65) 



L[h,J] = bxZ[h + iJ] - -J* U' 1 * J , (66) 

Here, we have introduced the notation Z[h] for the par- 
tition function of an ideal gas of molecules subjected to 
an applied field h, for which the total potential energy is 
Uo+U eK t[h]. In Eq. (|6"6"]) . Z[h + iJ] is the grand-canonical 
partition function for an ideal gas in which monomers of 
type i are subjected to a fluctuating complex field 



hi(v) ee hi(r) + iJi(r) 



(67) 



Hereafter, quantities such as Z and h that are defined in 
this ideal gas reference state will be denoted with a tilde. 

The functional derivative of lnZ[/i] with respect to h 
is the average monomer concentration 



(*(r)> 



8 In Z [h] 
Shi(r) 



(68) 



Here, (• • •) denotes an average taken in the ideal gas sub- 
jected to a field h. Higher derivatives yield 



(■») 



,r„) 



5^ \nZ[h] 



8h il (Y l )8h l2 (Y2) ■■■6h iri (r n ) 



(69) 
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where 

^W..;„( r l>"-> r n) =X] XI ( C arm 1 (ri)---C Qmi „(r„)) 
am si---s n 

(70) 

is an n-point intramolecular correlation function in 
the ideal gas reference state for monomers of types 
■ ■ ■ ,i n on the same molecule. 

A. Gaussian Fluctuations 



where Z s = e L ^ h > J 1 is the saddle-point value, and Zq is a 
factor arising from Gaussian fluctuations. The Gaussian 
contribution to the free energy in a homogeneous liquid 
is 

-lnZ G = \v [ lndet[G- 1 (k)U(k)] 

= \v { lndet[I + 0(k)U(k)] , (76) 
2 Jk 

where det[- • •] denotes the determinant of a 2 x 2 matrix. 



As already noted, a very simple form of mean-field the- 
ory is obtained by applying a saddle-point approximation 
to functional integral (|65[) . Requiring that 



6L[h, J] 
6 Mr) 



= 



(71) 



and using Eq. (|68|) yields the saddle-point condition 
given in Eq. (|40|) . in which /i,(r) denotes the saddle- 
point value of the field h = h+iJ defined in Eq. (|67|) , and 
in which (cj(r)} is the corresponding ideal-gas monomer 
concentration. 

A Gaussian approximation for Z is obtained by the 
approximating the deviation SL of L from its saddle point 
value by an expansion to second order in the deviation 
SJi(r) = Ji(r) — J/(r) of J from the saddle-point field 
J s . We thus approximate 



where 



(72) 



<V(ry) = - 6L[J ' h] 



5J t {r)dMr') 
= n i3 (v,v') + U^(v,v') , (73) 

or, in the Fourier space representation for a homogeneous 
fluid, 



fc 1 (k)=%(k)+?7 i 7 1 (k) 



(74) 



The "propagator" G obtained by inverting G _1 is the 
screened interaction potential identified by Edward s 11 ' 12 , 
which is closely analogous to the electrostatic screened 
interaction in the Debye-Huckel theory of electrolytes. 
It was noted by Edwards that results of the one-loop 
theory could be obtained from a perturbation theory in 
which monomers interact via this screened interaction. 
It was shown more systematically by one of usi£ that the 
diagrammatic perturbation theory that arises naturally 
from the auxiliary field method is equivalent to a type of 
molecular cluster expansion in which interaction lines (or 
"bonds") represent factors of the screened interaction. 
In the Gaussian approximation 



B. Incompressibilty and Regularization 

In the limit i; <C A, an explicit expression may be 
given for G(k) for a nearly incompressible homogeneous 
liquid. By taking the incompressible limit Bq — > oo , 
while approximating F(k/A) — 1 for k -C A, we obtain 



G«(k) 



i-2v X o\n(k)\n^(k) 

Cl + (k) -2v Xo Mk)\ 



(77) 



The T ' in the numerator denotes a contribution of unity 
to every matrix element. 

In the opposite high-wavenumber limit k 3> A, the 
bare potential U(k) in any liquid with a large but finite 
bare compression modulus Bq is assumed to become van- 
ishingly small, as a result of the decay of the crossover 
function F(k/A). This decay of the bare potential will 
also cause the screened potential G(k) to become very 
small for k 3> A. 

In Eq. ((76)) . the Gaussian contribution to the free en- 
ergy is expressed as a Fourier integral in which the inte- 
grand depends on G(k). This is a generic feature of the 
auxiliary field theory: Fluctuation corrections to mean- 
field results for all quantities of interest may be expressed 
as Fourier integrals involving factors of the screened in- 
teraction G(k). Because both £7(k) and G(k) vanish for 
k A, the use of any model with a nonzero range of in- 
teraction A -1 , rather than a point-like interaction, thus 
naturally introduces a cutoff length. 

The effects of a nonzero range of interaction A -1 may 
thus be crudely mimicked by treating A as a cutoff 
wavenumber, and simply suppressing contributions from 
wavevectors k 3> A in all Fourier integrals. This regular- 
ization scheme, which we will adopt, is equivalent to the 
use of a model in which U(k) is independent of k for all 
k < A, and zero for all k > A. In what follows, we will 
also restrict ourselves to the nearly- incompressible limit, 
in which G(k) is given for all k < A by Eq. (|77|) . 

We consider a Gaussian approximation for the free en- 
ergy density of a homogeneous liquid (e.i, either a binary 
homopolymer blend or a disordered diblock copolymer 
melt) in which the integral in Eq. (|T6[) is restricted to 
k < A, and in which the integrand is evaluated for k < A 
by considering the limit of large Bo. This yields a one- 
loop correction to the free energy density 



Z ~ Z S Z G 



(75) 



5f 



lnZ G /V 



(78) 
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pven by 



where 



2 



ln{[n + (q)-2« X o|n(q)|] ^} 



(79) 



|q|<A 



Here, f2+(q) and |f2(q)| are the sum of elements and de- 
terminant, respectively, of the ideal-gas correlation func- 
tion matrix fty(q), defined by analogy to Eq. I|35p for 
the related quantities fi+(q) and |0(q)j. 



V. FREE ENERGY IN BINARY BLENDS 

Several author a 5 ' 6 ' 7 ^ have obtained a UV divergent 
contribution to Xa(k = 0) in a binary homopolymcr 
blend by calculating the one-loop free energy Sf for a 
homogeneous mixture, as a function of composition, and 
then using Eq. (f3"8"| to extract x a (k = 0). In this section, 
we review and extend this approach. 



A. High-g Behavior of £1 and G 

To analyze the UV divergence of Sf, we will need an 
asymptotic expansion of the high-q behavior of the in- 
tramolecular two-point function Qjj(q), as an expansion 
in increasing powers of 1/q. In a binary blend of two 
Gaussian homopolymers, the function fiy (q) is a diago- 
nal matrix with elements 



fiy(q) = %CjA(q) 



with 



A(q) = N t g(Ql) 
where Q 2 = q 2 Nbf/6, and where 

g{x) = 2(e~ x - l + x)/x 2 



(80) 



(81) 



(82) 



is the Debye function. The required high-q expansion of 
Di (q) may be obtained by simply dropping the exponen- 
tially decaying term in the Debye function (which is not 
an analytic function of 1/q). This yields an approxima- 
tion 



A(q) 



A (0) (q) - ,, 



A (0) (q) 



2Ni 

Qi 



A (1) (q) 



12k 1 

V q 
_-72l 2 1 
v 2 Ni q 4 



(83) 

(84) 
(85) 



v/b 2 



where v is a monomer reference volume, and k 

A corresponding expansion of the screened interaction 
to the same order yields 



G<°>(q) = 

G (1) (q) = 



with I = 



^(q) 

Mi + <hh, 
z = 



v 

m 

V 



2V 7 



N 2 



and 



-(fell <l>i<hhh 
(frifohh -4>¥i 



(87) 



(88) 



In Eq. dHSJ and fEZ}, the quantities G(°)(q) and (?W(q) 
arc written with no monomer type indices i and j to 
indicate that the values of these quantities are actually 
independent of i and j . 



B. Free Energy Density 

In a binary homopolymer blend, Eq. (|79[) f° r ^/ re ~ 
duces to 

6f~ \ [H(n 1 + n 2 -2v Xo n 1 n 2 )B v] , (89) 



The dominant contribution to the argument of the log- 
arithm in Eq. (|89p in the high-fc limit arises from the 
leading order contribution to f2+ = ill + Q 2 , which is 



^ 0) (q)^E c ^ (0) (q) 



12/ 1 

v 2 q 2 



(90) 



Note that n { °\q) = 1/G< )(q). Factoring ^ B out of 
the argument of the logarithm yields an expression 



8f = 6fW + ± I In 



{l + G (0) [a« + -2vxo|^|]} (91) 



in which Af2 + (q) = £l + (k) — Cl^ (k) , and in which 
= i Jin (flpBov) 



12tt 2 



h '( 12So i) + 5 



A 



(92) 



is the most strongly UV divergent contribution. 

Upon expanding the integrand of the remaining inte- 
gral in Eq. (|9ip in powers of 1/q, we find that the leading 
order terms are of 1/q 2 , and yield UV divergent contri- 
butions of 0(A) to the integral, but that all subsequent 
terms in the expansion are UV convergent. The total UV 
divergent contribution to Sf is given by a sum 



Gij(q) — &°\q) + G^(q) + Grf'(q) 



(86) 



Sf ~ 5f {t)) + Sf {x> + Sf 



(x) 



f(i) 



(93) 
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where 



sf M = -xo G^\n^\ 



SfU = 1 j G (o)fjW 



(94) 
(95) 



where |^ (0) (q)l = c x c 2 D\ ' (q)L> 2 (q) and Q 
J2i c iD^\oi). Completing the integrals yields 

of (x> = ^--r^i^A 

IT V I 



n(i) 



(96) 



(i) 



3 ^f + M )A ( 97) 



2tt 2 vI V N x N 2 



The quantity 5f(°\ which is the one-loop contribution to 
the free energy density of a nearly incompressible blend 
of infinite chains with x = 0; arises solely from the strong 
repulsive interactions that suppress overall concentration 
fluctuations. This quantity generally depends on compo- 
sition, because of the composition dependence of /, and 
thus can contributes to Xa(0). It is independent of com- 
position only in the special case 1% = 1% = I. The con- 
tribution 5f^ x ' is a negative contribution that reflects 
a reduction in the average interaction energy from its 
mean-field value by correlations. We show in Sec. IXIII 
that 6/W arises from changes in packing near the chain 
ends. 



C. Direct Correlation Function 

Once the UV divergent contribution Sf is known, the 
corresponding contribution SCtj (k = 0) to the direct cor- 
relation function may be obtained from the relation 



sa 3 (o) = 



d 2 (Sf) 
dcidcj 



(98) 



Differentiation of Eqs. (|9"2] ) . ([55 ]) , and ([97 ]) yields a UV 
divergent contribution 



v 2 ld 



13 A 3 



^(0) 



6V klj ( (pill i 



2tt 2 I 3 V N l N 2 



A 



A 



(99) 
(100) 

(101) 



By applying Eq. (f3"6")> . we obtain a corresponding UV 
divergent contribution to Xa(0), 



Ol -^ A 3-SS XO A 



24tt 2 Z 2 



TT 2 P 



3 hhih-h) 



2tt 2 



P 



h_ 

Ni 



h_ 

N 2 



A ,(102) 



of the form given in Eq. (|41j) . Note that the coefficient 
A of the A 3 contribution and coefficients Hi and H 2 of 
the terms linear in A/Ni all vanish in the case h = l 2 of 
two polymers with equal statistical segment lengths. 

An alternative method of deriving £Gy(0), which is 
useful for comparison to the subsequent calculation of 
SCij(q) at q f 0, is to apply Eq. d32J) to Eq. 1(75]) 
before evaluating the Fourier integral. A straightforward 
differentiation yields an integral 



6C l3 (0) 



i d 2 



2 dcidcj Jci 



lndet[G _1 (q)U(q)] 



D i (q)G ij (q)G ii (q)£> i (q) . (103) 



We show in what follows that this expression can be re- 
covered by taking the k — » limit of our expression for 

JCy(k). 



VI. GAUSSIAN APPROXIMATION FOR S(k) 

In this section, we derive an Gaussian approximation 
for the inverse correlation function 5^ 1 (fc) at k ^ 0. The 
derivation is based on a calculation of the second func- 
tional derivatives of In Zq [h] with respect to an external 
field hi. 

For this purpose, we now adopt a compact notation 
for position and monomer type indices in which 5(1,2) 
is used as shorthand for a function Si x i 2 (ri, r 2 ). In this 
notation, an integer label / is used as shorthand for a 
position n and a monomer type index ij. For example, 
identity (fl8)) becomes 



5(1,2) 



S 2 \nZ 
5h(l)Sh(2) 



(104) 



where 6/Sh(2) denotes a functional derivative with re- 
spect to hi 2 (r 2 ). 

The one-loop contribution to 5(1,2) is given by a 
derivative 



55(1,2) 



5 2 lnZ G 



8h(l)5h(2) 



(105) 



in which Zq is the Gaussian contribution to the partition 
function of a system in an applied field h. The Gaussian 
contribution to the free energy of a system that is sub- 
jected to an inhomogeneous field h is a functional 



InZr 



1. 



-TrlnG -1 ?/ 



(106) 



Here, G _1 (l, 2) is a integral operator, which is generally 
not translationally invariant in a fluid that is subjected 
to an inhomogeneous external field h. The symbols In 
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and Tr denote the generalized logarithm and trace for 
such an operator. The function 

G~ l (l, 2; [h]) = & 2) (1, 2; [h]) + CT^l, 2) (107) 

is a functional of the field h, which is the saddle point 
of h + i J obtained in the presence of an external field h. 
The Gaussian integral Zq is thus also a functional of h. 

A straightforward differentiation of Eq. (|106[) twice 
with respect to the external field h yields 



cJS(l,2) = 



S 2 \nZ G [h] 5h(l')5h{2') 
v jv 5h{V)5h(2') Sh{l) 5h{2) 
SlnZ G [h] Sh(l') 



n- Sh(V) Sh(l)6h(2) 
Here, we have introduced the shorthand 



(108) 



(109) 



for integration over a position rj and summation over 
allowed values of a corresponding monomer type index 
ij. To evaluate Eq. (|108|) . we must evaluate both the 
first and second derivatives of Zq [h] with respect to the 
saddle point field h, and the first and second derivatives 
of h with respect to the external field h. 

To evaluate the required derivative of Zq [h] , we apply 
the identities 13 



^TrlnG" 1 
5h(l) 



SG~ 1 (2,3) 



G(2,3)- 
i 8h{l) 



^(1,-,") = ^)(l,...,n,n + l) (110) 
6h{n + l) v.j. J ) 



6G{1,2) 
5h(3) 

to obtain 

5\hZq 



G(1,1')^ (3) (1'.2',3)G(2',2) 



1' J 2' 



Sh(l) 



ft (3) (l,2,3)G(2,3) (111) 



2 J3 



and 



5 2 lnZ G 
5h{\)8h{2) 



f* (4) (l,2,3,4)G(3,4) 



(112) 



3 Ji 



ft (3) (l, 3, 4)G(3, 5)0 (3) (2, 5, 6)G(6, 4) 



3 J4 J5 J6 



To evaluate the derivatives of h with respect to the 
external field h, we differentiate the saddle saddle-point 
condition 

h(l) = h(l)-J U(l,2)m (2) , (113) 



with respect to h. Note that fl^(2) is a functional of 
h(l), with derivatives given by Eq. (|110p . A single func- 
tional derivative yields 



Sh(2) 
2 Sh(l) 



5(2,4)+ / fi( 2 >(2,3)[/(3,4) 



5(1,4) (114) 



Let 



5- 1 (l,2) = ^ 1 (l,2) + [/(l,2) 



(115) 



denote the inverse correlation function obtained in the 
mean-field approximation, in which G(l,2) = — [7(1,2). 
Using this definition for 5 _1 and Eq. ([75)) for G^ 1 , we 
may rewrite Eq. (|1 14[) in either of the equivalent forms 



Sh(2) 

dh(T) 



5(1, 3)^(3, 2) 
U-\l, 3)G(3, 2) 



(116) 
(117) 



A second functional derivative yields 

6 2 h(3) [ [ [ [ ^M 1 ') Sh(2') 

8h(l)6h(2) Z J v J v J 3 ,] v 5h{l) 5h{2) 

x^ 3 )(l',2',4')C/(4',3')^(H8) 

An explicit expression for SS (1,2) may then be ob- 
tained by combining Eqs. (fTOg]) . (fTTT]) . (fTTBl) or ([TT7|l . 
and pi8[) . We are more interested, however, in obtaining 
a one-loop contribution to the inverse correlation func- 
tion S" -1 (1,2) that appears in the OZ relation. To the 
order required here, this is related to 55(1, 2) by 



55^(1,4) = - / / 5(1,2)55(2,3)5(3,4) (119) 

J 2 J 3 



Combining Eq. (| 1 1 0|) with our expression for 55(1,2) 
yields a one-loop contribution 



55- x (l,2) 

in which 

5A(1,2) 



n- 1 (l,l')5A(l',2')n- 1 (2',2) 



5 2 In Z G 
Sh(l)Sh(2) 

^ 3 )(1,2,3)G(3,4)^ 
.; ./ 1 5/i(4) 



(120) 
(121) 



The required functional derivatives of In Zq are given in 
Eqs. (fTTTj) and fTT2")) . 



VII. DIAGRAMMATIC APPROACH 

It is convenient to introduce a diagrammatic repre- 
sentation of cluster integrals such as those obtained in 
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(a) 



(b) 



(c) 



FIG. 1: The three one-loop diagrams of Q vertices and — G 
bonds that contribute to the expression for SA(1, 2) given in 
Eq. (|12ip . Diagrams (a) and (b) represent cluster integrals 
arising from the first and second terms of the r.h.s. of Eq. 
(fTT2l for 5 2 lnZ/5h(l)8h(2), respectively, in the first line of 
Eq. (|12ip . Diagram (c) represents the second line of Eq. 
(fl2Tl) , in which S\nZ/5h{4) is given by Eq. (TTTTj) . 



Eqs. (fTTTj) . (fTT2)) . and (fmj) . We adopt the diagram- 
mar used previously in Ref. [13]. Some examples of the 
type of diagrams used here are shown in Fig. [1] Each 
integral is represented as a diagram containing vertices 
and bonds. An n-point vertex, shown as a shaded circle 
with n smaller circles around its circumference, repre- 
sents a function of n coordinates and type indices, such 
as f2' n )(l, . . . , n). Each of the small circles around the 
perimeter of a vertex is either a field circle (shown black- 
ened), which represents an argument of the correspond- 
ing function that must be integrated over, or a root circle 
(white), which represents a fixed parameter, rather than 
a integration variable. Each bond represents a function of 
two coordinates and type indices, such as the bare inter- 
action {7(1,2) or the screened interaction G(l,2). Each 
bond must be connected at each end to either a vertex 
field site, or a free root site. (A free root site is small 
white circle that is not associated with a vertex, which 
is used simply to indicate the arguments associated with 
the free end of a bond are known parameters). The value 
of a diagram is the value of the integral obtained by in- 
tegrating over the coordinates associated with all of the 
black circles, divided by a combinatorical prefactor that 
is given by the order of the group of permutational sym- 
metries of the diagram. The diagrams discussed here are 
all diagrams of Q- vertices and — G bonds, in which a fac- 
tor of Qy 1 ' is associated with each n-point vertex and a 
factor of — G is associated with each bond. In diagrams 
with —G bonds, vertices with no white circles (i.e., with 
no root circles) must each have three or more black (field) 
circles, representing factors of with n > 3. 

Fig. [T] shows a diagrammatic representation of Eq. 
(fT2l"j) . for 5A(1,2), in which we have used Eqs. (fTTT]) 
and (|112[) for the required functional derivatives of In Zq- 
The correspondence between diagrams and integral ex- 
pressions is discussed in the figure caption. 

It was shown in Ref. [13] how In Z and various correla- 
tion functions could be expressed to any order in pertur- 
bation theory as sums of well-defined infinite sets of such 



cluster diagrams. It was found that the two-point cor- 
relation function S(l,2) can be expressed as an infinite 
sum 



5(1,2) 



Sum of all connected diagrams 
of O vertices and — G bonds 
with 2 roots circles labelled 
1 and 2 



(122) 



The intramolecular correlation function 0(1, 2) in an in- 
teracting fluid is given by the subset of these diagrams 
in which both of the root sites are one the same vertex 
(the root vertex): 



0(1,2) = 



Sum of connected diagrams 
of O vertices and — G bonds 
with 2 roots circles labelled 
1 and 2 on the same vertex 



(123) 



Among the diagrams described in Eq. (|123[) is one con- 
sisting of a single two-point O vertex, with no bonds, 
which represents the ideal gas contribution 0(1,2). The 
corresponding expansion of the intramolecular correla- 
tion function o (l, 2) for molecule of a specific type in a 
mixture may be obtained by replacing the O root vertex 
in each diagram of Eq. (|123[) (i.e., the vertex with two 
white circles) by a corresponding O a vertex, representing 
a factor of the intramolecular correlation function O a for 
the specified molecule type, while using SI vertices for all 
other vertices in the diagram. 



A. Diagrammatic Resummations 

The screened interaction — G used in the cluster ex- 
pansions described above can be expressed algebraically 
as an infinite geometrical series 



g = -\i + u*n] 



* u 



-u + u *n*u - 



(124) 



where O denotes &( 2 \ This can also be expressed dia- 
grammatically as the sum of an infinite series of all pos- 
sible chain diagrams of alternating — U bonds and fl^ 
vertices. By substituting this diagrammatic expansion of 
— G into the above expansions of S and O, we may ob- 
tain a formally equivalent expansions of these quantities 
in terms of diagrams of O vertices and —U bonds. That 
is, the perturbation theory may be expressed in terms of 
either the screened interaction G or the underlying bare 
interaction U. The descriptions of the infinite sums of 
diagrams of O vertices and —U bonds required to con- 
struct S and O are identical to those given in Eqs. Q122[) 
and (|123p for diagrams of O vertices and — G bonds, ex- 
cept for a replacement of — G bonds by — U bonds, and 
a change in the rule for the nature of the allowed O ver- 
tices: Two-point f)( 2 ) vertices with two field circles and 
no root circles are allowed in diagrams of — U bonds, but 
are prohibited in diagrams of — G bonds. 
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A Dyson equation for 5(1, 2) may be obtained by defin- 
ing a function 

{Sum of bond irreducible diagrams *| 
of vertices and — G bonds with > 
2 roots circles labelled 1 and 2. J 

(125) 

A "bond irreducible" diagram is one that cannot be di- 
vided into two disconnected pieces that each contain one 
of the two root circles by cutting or removing only one 
bond. Eq. (|125p can also be expressed as a sum of all 
bond-irreducible diagrams of O diagrams and — U bonds, 
if fj( 2 ) vertices with two field sites and no root sites are 
allowed. This set of diagrams includes the trivial diagram 
consisting only of a 0( 2 ) vertex with two root sites and no 
bonds, which represents the function 0( 2 )(1, 2). Thus, to 
a first approximation (or "tree level"), A(l,2) ~ 0(1,2). 
The only one-loop diagrams that contribute to A(l, 2) are 
the three diagrams shown in Fig. [TJ The expression for 
SA(1, 2) given explicitly in Eq. (|12ip is thus the one-loop 
contribution to the quantity A(l, 2) defined in Eq. (|125p . 

The function S can be expressed in terms of A as a 
geometric series 

S = A-A*U*A + --- 

= A-A*U*S , (126) 

or diagrammatically as a sum of chain diagrams of A 
vertices connected by — U bonds. Resumming this series, 
or solving the recursion relation, yields 

S _1 (l,2) = A- 1 (l,2) + [/(l,2) (127) 

Approximating A(l, 2) by 0( 2 )(1, 2) yields the mean field 
approximation for S. Note that Eq. (fT2T)l for ,S _1 (1,2) 
is not the same as the generalized OZ equation. 

B. One-Loop Contributions 

To calculate one-loop corrections to S^ 1 , it is useful to 
define 

A(1,2) = 0(1,2) + <5A(1,2) , (128) 

where <5A(1,2) is a sum of all contributions to A(l,2) 
other than the tree-level contribution 0(1,2). To first 
order in a loop expansion 5A(1, 2) is given by the sum of 
the three one-loop diagrams shown in Fig. [1] Substitut- 
ing Eq. (I128|) into Eq. (|127|) yields a geometrical series 
in which the first few terms are 

s- 1 = u + n- 1 - fr 1 *5A* fr 1 + • •• (129) 

The one-loop contribution to S'~ 1 (l, 2) is thus given by 

ss- 1 = -fr 1 * sa * fr 1 (130) 

with 5 A approximated by the sum of the three one-loop 
diagrams for A. 



The one-loop contribution to the OZ expression for 
may also be expressed as a sum 

SS' 1 = 5Q- 1 - SC (131) 

where <50 _1 and SC represent one loop contributions to 
_1 and C, respectively. The one-loop correction to _1 
is given by the convolution 

sn- 1 ~ -fr 1 * sn * fr 1 (132) 

where S£l represents the one-loop correction to O. Ac- 
cording to the diagrammatic rule given in Eq. (I123|) , the 
one-loop contribution to 0(1,2) is given by the sum of 
diagrams (a) and (c) in Fig. Q] It follows that we may 
identify — SC with the remaining contribution to Eq. 11301 
that arises from the contribution of diagram (b) to SA. 
That is, to first order in a loop expansion, 

sc ~ fr 1 * s * rr 1 , (133) 

where S is the value of diagram (b). This expression was 
obtained previously^ as the one-loop contribution to a 
general diagrammatic expansion of the direct correlation 
function. 



1. Intramolecular Correlations 
The quantity <50y(k) may be expressed as a sum 

$n -(k) = $^<yfi o ,«(k) (134) 

a 

where Sfl a ,ij (k) is a one loop correction to the intramolec- 
ular correlation function Q a ,ij (k) for molecules of species 
a. Using the diagrammatic rules discussed above, the 
one-loop contribution Srt a ^j (k) may be expressed as a 
Fourier integral 

SCl a>ij (k) = p a [ I a % (k) - u> a>ij (k)G ab p b I D 0) } (135) 
Here. 

Ok) = -^<^(k,-k,q,-q)G«(q)(136) 

where u) a ,ijW = ^a,i( k > -k )' anci 

Gab = N ^Gu (k = 0)JV« (137) 

kl 

where A a j. is the number of monomers of type A: on a 
molecule of type a. Here and in what follows, we use the 
notation 




(138) 
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for Fourier integrals. To obtain the above expression, 
we have used the identity (D[ 3 ^- fe (k, — k, 0) = NakCo^ (k). 
The first term in square brackets in Eq. (|135[) corre- 
sponds to the diagram (a) in[IJ while the second corre- 
sponds to diagram (c). To obtain the above expressions 
for 6il at ij(k), the root vertices in these diagrams must be 
taken to be fl a vertices. A corresponding contribution for 
the sum <5f2jj(k) may be obtained either by adding the 
results for <5f2 0j jj for different species in a mixture, or by 
taking the root vertices in these diagrams to be fl vertices 
(with no species index). 



2. Direct Correlation Function 



Eq. (|133p for 8C may be expressed more explicitly in 
Fourier space as a product 



sCijiv-n-^k^Mrij. 1 ^) , (139) 



in which 



£y(k) 



n^(k,q_,-q + )G w (q+) (140) 



x ng(-k,-q_,q + )<3« n (q_) , 



is the Fourier transform of diagram (b) of Fig. 1 , where 
q± = q± |. 



VIII. CANONICAL VS. GRAND-CANONICAL 
ENSEMBLE 

Because we have worked thus far in grand-canonical en- 
semble, the one loop contribution to S^ 1 derived above 
is a correction to the mean-field result for a system with 
a fixed set of chemical potentials, rather than fixed set 
of molecular concentrations. In grand-canonical ensem- 
ble, a one-loop approximation for the free energy will 
generally yield a slightly different concentration for each 
type of molecule than that obtained at the same chem- 
ical potentials from the saddle-point approximation. In 
this section, we first calculate the difference between the 
molecular concentrations obtained in the mean-field and 
one-loop approximation at equal chemical potentials, and 
then use this to obtain a one- loop approximation for S^ 1 
for a system with fixed molecular concentrations. 



A. Concentration at Fixed Chemical Potential 

In grand-canonical ensemble, molecular concentration 
is given by a derivative 



Pa 



1 d\nZ 
V dp a 



(141) 



The concentration obtained in a Gaussian or one-loop 
approximation thus differs from that obtained in mean- 
field theory at equal chemical potential by an amount 



5p a 



1 dhxZ G 

V dfl a 



(142) 



Here, the derivative is evaluated at fixed temperautre and 
fixed values of the chemical potentials of species other 
than a. Because the saddle-point field h will shift in 
response to changes in chemical potential, however, and 
In Z G is expressed as a functional of the saddle-point field, 
this derivative is a sum 



VS Pa = 



d\nZ G 



0p a 



dh{l)8\nZ G 



1 1 dfi a 5h(l) 

where Sin Z G /Sh(l) is given in Eq. (jllip . and where 
d\nZ G 



(143) 



dPa 



n a (l,2)G(2,l) 



(144) 



1 J2 



A straightforward functional derivative yields the iden- 
tity 



dh(l) 

dfl a 



G(l,2)OM(2) 



(145) 



where Q^) (r) is the contribution of molecules of type a to 
the concentration of monomers of type i in the reference 
ideal gas. This is given in a homogeneous fluid by a 
constant 



(146) 



for all r, where N a i is the number of monomers of type 
i on a molecule of type a. By combining Eqs. (I143|) and 
(|144j) . we obtain 



V6 Pa 



1 

2 ji Ji 
1 

2 



Q (1,2)G(2,1) (147) 
^(1)6(1, 2M2,3,4)G(3,4) 



1 J2 J3 J A 



This expression is shown diagrammatically in Fig. [2] 
Evaluating the Fourier representation of these diagrams 
yields 



(148) 



where I a °^ and G a b are defined in Eqs ()136|) and (|137|) . 



B. S(k) at Fixed Concentration 

We would like to obtain an expression for the one-loop 
correction to S(k) in an interacting liquid with a fixed 
concentration p a for each species of molecules. To do 
so, we consider a one-loop approximation in which each 
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(a) (b) 



FIG. 2: The two one-loop diagrams of £2 vertices and — G 
bonds that contribute to the molecular number density p a - 
Diagrams (a) and (b) correspond to the first and second lines 
of Eq. (|147jl , respectively. Here, whitened vertices (larger cir- 
cles) are used to represent factors of the ideal-gas intramolec- 
ular correlation function Q a for a specific molecular species a, 
while the gray vertex in diagram (b) represents a factor of fi, 
as defined in Eq. ((22)) . 

chemical potential is shifted by an amount Sp a from the 
value used in the mean- field calculation, where 5p, a is 
chosen so as to produce a shift — Sp a that exactly cancels 
the the one-loop contributions to p a at fixed chemical 
potential. To first order in a loop expansion, the only 
effect of this shift in chemical potentials is to shift the 
intramolecular correlation function by an amount 

5n aiij (k) = -Sp a Q atij (k) (149) 

where Sp a is the one-loop contribution to p a at fixed 
chemical potential. There is no corresponding correc- 
tion to CV/(k), because this is given to zeroth order by a 
quantity — £/jj(k) that is independent of p. 

Upon combining Eq. (|149|) with Eq. (|148[) for 5p a , we 
find that the contribution of diagram (c) of Fig. [T]to 5Vt a 
is precisely cancelled by the contribution of diagram (b) 
of Fig. [5] to 5 p a - The contributions arising from diagram 
(a) of Fig. Q]and diagram (a) of Fig. [5] yield a correction 
Stt^ijOi) = p a 5u) aj ij(}s), where 

Su^a (k) = /£> (k) - (k)/(°) (iso) 

This one-loop contribution to Lo a ,ij(k) at fixed molecular 
concentration may also be expressed as an integral 

^S( k )^-jE / ^«(k,-k,cb-q)G«(q) , 
ij Jq 

(151) 

where 

&(k,-k,q,-q) = <^(k, -k, q, -q) 

-43i(k)<L(q) (152) 

(2) 

This expression for SutfQn) was obtained in Ref. [13] by 
similar reasoning. It was obtained previously by Barrat 
and Fredrickson by calculating a one-loop correction to 
LUij(k) in canonical, rather than grand-canonical, ensem- 
ble. 



IX. ASYMPTOTIC ANALYSIS OF 
SINGLE-CHAIN CORRELATIONS 

To show that the one- loop approximation for ^"^(k) 
is renormalizable, we must show that the UV divergent 
parts of the one-loop contributions to uiij (k) and dj (k) 
each have a specific wavenumber dependence. In both 
cases, the proof relies critically upon an analysis of the 
dominant asymptotic high-q behavior of a 3- or 4-point 
intramolecular correlation function. To show that domi- 
nant UV divergent part of Eq. (|139[) for SCij (k) is inde- 
pendent of k, as required by our criteria for renormaliz- 
ability, we must examine the asymptotic high-q behavior 
of the 3-point function f2^(k, q , — q + ) that appears in 
Eq (|140[) . Similarly, to show that the UV divergent part 
of Eq. (|15ip for (5aj a .r/(k) has the non-trivial wavenum- 
ber dependence implied by Eq. (|48p . we must exam- 
ine the asymptotic behavior of the four-point function 
-0^ ; (k, — k, q, — q). In both cases, we need the asymp- 
totic behavior of a multi-point correlation function in the 
limit of large integration wavevector q, for arbitrary k. 
The correlation functions required here are closely related 
to, but distinct from, the 3- and 4-point vertex functions 
introduced by Leibler-21, which are defined as the 3rd 
and 4th functional derivatives of the SCFT free energy 
-F[(c)]. For simplicity, we limit ourselves in this section 
to the identification of the dominant high-q behavior for 
homopolymers. A more systematic asymptotic expan- 
sion for homopolymers and the generalization to diblock 
copolymers are given in the appendix. 



A. Multi-Point Correlations 

To begin, we consider some general features 
of the rt-point intramolecular correlation function 
d)( n )(ki, . . . ,k n ) for a Gaussian homopolymer. This 
quantity is given by an integral 

N 

cD(")(k l5 ...,k n ) = fd n s^e^i^ , (153) 
o 

where (• • •) indicates an average over conformations of 
a single Gaussian chain. Here, we have introduced the 
notation 

N N N N 

J d n s = J ds n - J ds 2 J dsi (154) 



for an integral over n contour variables. Because we con- 
sider only homopolymers in this section, no monomer 
species indices are needed, or used. 

To calculate a)'™'({k}), we divide the integral over 
< si , S2, ■ ■ ■ , s n < N into n! contributions arising 
from different ways of ordering the values of s\, S2, 
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. . . , s n , each corresponding to a different ways of or- 
dering labelled monomers along a chain. We first con- 
sider a restricted n dimensional integral in which the 
variables si, . . . , s n are restricted to the ordered subpace 
S\ < S2 < • • • < s n , and then reconstruct the original in- 
tegral by summing over all n\ permutations of the index 
labels i — 1, . . . , n. By this method, Eq. (|153[) may be 
rewritten as a sum 



^" ) (k 1 ,...,k„) = ^^(kf,...,k^ 

{P} 



(155) 



in which the list kf , . . . , k^ is a permutation P of the 
original list of wavevectors ki, . . . , k„, and J2{p} denotes 
a sum over all possible permutations of the n wavevec- 
tors. Here, 



JY 



w< n >(k 1> ...,kn) = J d n s (e'^VR^)) (i5 6 ) 



is an ordered integral, in which we introduce the notation 

N N s 3 s 2 

J d n s = J ds n ■■■ J ds 2 Jd Sl , (157) 



to indicate an ordered integral over the space < si < 

S 2 < S 3 < • • ■ < S n < N. 

To obtain an expression for the ordered integral , 
we first rewrite the sum of dot products that appears in 
the exponential in Eq. (|156p as a sum 

n n— 1 

J2 k j • R ( s j) = e « • H(«n) " e ' AR J+U ( 158 ) 



in which 



and in which 



ARy = R(s,) - R( Sj ) 



(159) 



(160) 



where ki , . . . , k„ denotes the ordered list of arguments of 
the ordered integral uj^ (ki , . . . , k n ) . Using the property 



^e-[R0)-R(s')]\ _ -e 2 b 2 \s-s'\/6 



(161) 



of a Gaussian chain, where e 2 = |e| 2 , we then obtain an 
explicit expression 



^ n) (k!,...,k n )= J Ae-^i'l^ 2 ' 2 ^^^ (162) 



where 



is a difference in monomer contour variables, for any set 
of wavevectors ki, . . . ,k„ for which e„ = ^™=i k i = 0> 
as required by translational invariance. 

By changing variables in the above integrals to contour 
variables = Si/N, for which < I, < 1, we find that 



(k u ...,k n ;N,b) = JV n wW(K 1 ,...,K n ) 



^ n \k u ...,k n ;N,b) = iV"^ ( " ) (K 1 ,...,K„)(164) 

where the functions uj^ and ur- n ' depend only upon the 
re-scaled wavevectors K.,- = kj^/Nb/6, and 



(165) 



&r(Ki,...,K n )= / (fse 



Here, £" 2 = |Ki + • • • + K,-| 2 , and the integral in Eq. 
(|165p is taken over < Si < §2 < ■ ■ ■ < s n < 1. The 
function u>( n ' is an n-point generalization of the Debye 
function, which is related to u/^ by a sum over permu- 
tations analogous to Eq. (|155l) 



B. Three Point Function 

Consider the asymptotic high-g behavior of the func- 
tion w( 3 )(K, Q , — Q+), which is needed in Eq. (|140j) to 
calculate 8Cij{K). Adding the 3! = 6 permutations of 
the monomers or wavevector labels yields 

^ 3 )(K,Q_,-Q+) = 2^ 3 >(K,Q_,-Q + ) 
+ 2^ 3 >(K,-Q+,Q_) 
+ 2w (3) (Q_,K,-Q+) (166) 

where Q± = Q ± K/2. Equal contributions to this sum 
are made by permutations that are related by a reversal 
of the order of the wavevector arguments, so that, e.g., 
^ 3 >(K,Q_,-Q + )=^ 3 )(-Q+,Q_,K). 

We are interested here in the asymptotic behavior of 
&( 3 )(K,Q_,-Q + ) in the limit Q 2 > 1 and Q 2 > K 2 , 
for otherwise arbitrary K. First, consider the first two 
ordered integrals in Eq. (|166[) . in which K is the first 
of the three arguments. These may be related by taking 
Q — > — Q, so we need consider only <j^(K, Q_, — Q+). 
This is given by an integral 



1 

^ 3 >(K,Q_,-Q + ) = J d 3 se- Q + s ™- K2s ^ 



(167) 



S%j — S{ Sj 



(163) 



In the limit Q 2 3> 1 of interest, this integral is dominated 
by contributions in which the separation S32 is of order 
l/Q 2 , yielding S32 <C 1. The dominant asymptotic be- 
havior in this limit may thus be obtained by replacing the 
integral with respect to S3 over the domain §2 < s 3 < 1 
by an integral over the seminfinite domain §2 < S3 < 00, 
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while still requiring that 0<si<S2<l. To leading or- 
der in an expansion in powers of l/<5, this approximation 
yields 



to 



(3) 



(K,Q _,-Q+) 



where 



g(K 2 ) = 2 / ds 2 / dsi e 



(168) 



(169) 



is the Debye function. Because the result is in- 
variant under Q — * — Q, this approximation yields 



•,(3) 



(K,-Q + ,Q_) 



,(3) 



(K, Q-,— Q+). Next, we con- 



sider the remaining ordered integral, 



,(3), 



Q_,K,-Q H 



d 3 s e -Q 2 + s ™-Q- s ^ 



(170) 



By reasoning similar to that discussed above, we may 
approximate this integral in the limit Q\ > 1 by an 
integral over a domain — oo < si < ,§2, S2 < S3 < oo, and 
< §2 < 1. This yields a leading order approximation 

<^ 3 )(K,Q_,-Q + )~ i J + G(Q- 6 ) , (171) 

which does not contribute to the leading order 0(Q~ 2 ) 
term in ui^ 3 \ The leading order contribution to ili^ thus 
arises from the four permutations in which K is either 
the first or last argument of cD^ 3 -* , and is equal to 4 times 
the r.h.s. of Eq. l(Igg]) . 

A more systematic asymptotic expansion, which is out- 
lined in the appendix, yields 



;( 3 )(K,Q_,-Q H 



£,(3,0) + fl(3,l) 



with 



',(3,0) ^_ 



,(3,1) _ 



^-L[^ 2 (3 + 4 a 2 ) 3 (X 2 )-4] 



(172) 

(173) 
(174) 



where a = Q-K./QK . Further terms in the expansion are 
not needed to analyze the UV divergence of 5Cy (k) . The 
observation that the leading order contribution a/ 3,0 ) has 
the same -ftf-dependence as the two-point function Cj^ — 
g(K 2 ) plays an essential role in our analysis of <5CV, (k). 



C. Four Point Function 

The function ip^ = tp^ (— k, k, — q, q) may be ex- 
pressed as an integral 



JV 



V> (4) = / d n s (( e lk - AR21 e »i- AR «) - ( e ik.AR 21 ^ e iq.AR 43 ^ 

(175) 



in which the integral ranges over < Sj < N for all i = 
1, ... ,4. Like the integral for this integral may be 
divided into contributions arising from each of the 4! = 24 
different possible permutations of the contour variables 
§1, . . . , S4. We thus define functions ip analogous to the 
ordered integrals ujj 4 ^ , in which arguments ±k and ±q 
are listed in the same order as order of the values of the 
associated contour variables. For example, 

^ 4 >(-k,-q,q,k) 

N 

= J jn g ^gik-AR.41 giq-AR.32^ _ ^gtk-AR 4 i ^giq-AR 32 ^ 


where the integration is over a subspace < s\ < S2 < 
S3 < S4 < N. We also define functions 

^ 4 )(K,-K,Q,-Q) = 7V-\; (4) (k,-k,q,-q) 

^ (4) (Ki,K 2 ,K 3 ,K 4 ) = iV-V 4) (k 1 ,k 2 ,k 3 ,k4)(176) 

analogous to Cj and Q_, which depend only upon the 
dimensionless wavevectors K = k-v/ Nb 2 /6 and Q = 
q^Nb 2 /6. 

"(4) 

The function ip is unchanged by permutations of its 
arguments that reverse the order (i.e., Ki,...,K4 — > 
K4, ...,Ki) or that reverse all of their signs (i.e., 
Ki, . .. ,K4 — > — Ki, . . . , — K4). Using these symmetries, 
we find that 

</> (4) = 4£(K, Q, -K, -Q) + 40(K, -Q, -K, Q) 
+ 2^(K,Q,-Q,-K) + 2^(K,-Q,Q,-K) 
+ 2±(Q, K, -K, -Q) + 2l(-Q, K, -K, QJ177) 

"(4) 

The values of %p arising from the 8 permutations in 
which ±K are the first two or the last two arguments, 
such as ip^ (K, — K, Q, — Q), vanish because of a cancel- 
lation of the first and second terms in the integrand of 
Eq. ((175]) . 

In the limit Q 2 > 1 and Q 2 > K 2 of interest, the 
dominant contributions to arise from the 4 permu- 
tations in which ±K are the 1st and last arguments of 
ijj. (The analysis required to show this is given in the 
appendix) . The sum of these permutations yields a dom- 
inant contribution of 0(Q~ 4 ) in an expansion in powers 

-(4) 

of 1/Q. For example, the function if) (k, q, — q, — k) is 
given by an integral 



<A (4) (K,Q,-Q,-K) 



d 4 s e 



43— W+ 



S32—K I21 



(178) 



Throughout this subsection, Q\ 



|Q±K| 2 . In the limit 



Q- ^> 1 of interest, the dominant behavior of this integral 
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^ (4) (K,Q,-Q,-K) 



-2Ka ±{Ka) 



dh / dh §4,1 e~ S41 (179) 



where a = K • Q/QK. The integral in the second line 
may be expressed as a derivative of Eq. p69[) for the 
Debye function, since 



dg(K 2 ) 

d{K 2 ) 



1 s 2 

= -2 J d§ 2 J dh hi e- R2s21 



(180) 



Upon adding the 4 permutations in which ±K are the 

- (4) 

first and last arguments of ip , the terms linear in a 
cancel. This yields a leading order contribution 

The dominant contribution to the function i// 4 ' = 
iV 4 ^( 4 ), may be expressed in terms of a partial derivative 



^(k,-k,q,-q) 



288a 2 d^ 2 >(k) 
~^U 2 dWV 



(182) 



of d>( 2 ) respect to b 2 . Here, we have introduced the nota- 
tion *!/)( 4 ' ) to denote the leading order term in an expan- 
sion of i/)*- 4 -* in powers of 1/q. This asymtotic expression 
is used in Sec. [X]to relate the one-loop correction to uj^ 
to a renormalization of the statistical segment length. 



G(q) oc q 2 and that Eq. 
large q. We thus expect 



!p2]) yields ^ cx <T 4 for 



d 3 q 4 i] Gi 



A 



(184) 



<3<A 



All subdominant terms in an expansion of the integrand 
in powers of 1/q 2 yield contributions to the integrand 
that are smaller by by factors of 0(l/q 2 ), which lead to 
UV convergent contributions to the integral. 

The ultraviolet divergent part of the one-loop approx- 

(2) 

imation for ou)\ is thus given by a cutoff integral 



&,,'*> (k) 



(4,0) 



(k,-k,q,-q)G(°)(q) (185) 



in which f// 4 ' ) denotes the leading order approximation 
for given in Eq. (TIE?)) , and G (0) (q) denotes the lead- 
ing order asymptotic approximations to Gjj (q) given in 
Eq. (|87|) (which is the same for all i and j). A straight- 
forward integration yields 



where 



<M 2) (k) 



d(b 2 ) 



5(b 2 ) 



2/ 2 



(186) 



(187) 



is a shift in the value of bf. 

The UV convergent correction to w(k) is given by the 
difference 



^.o) G (o) 



(188) 



X. CORRELATIONS IN BINARY BLENDS 

We can now analyze the UV divergent contributions 
to S^j (k) in a binary homopolymer blend. We consider 
the intramolecular and direct correlation functions sepa- 
rately. 



A. Intramolecular Correlations 

The one-loop contribution to cu a y (k) in cither a blend 
or a copolymer melt is given by Eq. (|151[) . In case of a 
binary homopolymer blend this simplifies to 



fcf } (k) 



^ (4) (k,-k,q,-q)G ri (q) (183) 



between the one-loop integral expression and the UV di- 
vergent part. Here, ip^ and -0( 4 '°) represent the exact 
and leading asymptotic expression for ?/>( 4 ) (k, — k, q, — q), 
respectively, while Gu and G^ represent the exact and 
leading asymptotic expressions for G M (q) . This is a con- 
vergent integral, which we evaluate by taking A — > oo. 
It is straightforward to show that both terms in this in- 
tegrand may be expressed as vN 3 (where N is Ni or 
7V2) times a function of a dimensionless wavenumber qR 
(where R = y/Nb and b is b\ or 62) and of the dimen- 
sionless parameters XoV, N1/N2, and 61/62- By non- 
dimensionalizing the measure in the wavevector integral 
by length scale R, we may express Sm* as a convergent 
dimensionless integral times a prefactor of vN 3 /R 3 = 
N 2 /N 1 ' 2 . The corresponding non-dimensional form of 
Sfl* is given in Eq. (|5i"j) . 



where the integral is constrained to |q| < A, and where 
we have introduced the notation uj^ = oj^ and ipi^ = 
^itii appropriate for a homopolymer mixture. 

To determine the order of the dominant ultraviolet 

(2,) _____ 

divergence of 8u>\ ; (k), we note that Eq. (|57)l yields 



B. Direct Correlation Function 

The one-loop contribution to the direct correlation 
function is given by Eqs. (|139|) and (|140|) . In case of a bi- 
nary homopolymer blend, for which all of the monomer 
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type indices of the functions Q$ or Cl^ft must be the 
same, Eqs. (|139j) and (|140p may be combined and sim- 
plified to obtain 



K7y(k) 



where 



Here, fL 



i /" D i (k,q)Gy(q+)Z? i (-k,-q)G <i (q_) , 

(189) 

(190) 



A(k,q) = 



_ ^ 3) (k,q_,-q+) 



fil 2) and n (3) 



)(3) 



f2^ are the intramolecular 
two- and three-point functions for homopolymers of type 
i, respectively. 

Consider the k — > limit of Eq. (|189p . By using the 
long-wavelength limits 



n 

^ 3) (0,q, 



(2) 



(0) = A.c, 
-q) 



we find that 



lim Z?i(k, q) 



f2i(q) 



(191) 



(192) 



By substituting this into Eq. (|189[) . we may immediately 
confirm that k — ► limit of Eq. (flM)) for <5Cy(k) is 
equivalent to Eq. (| 1 03[) for (5Gy(k = 0), which was ob- 
tained by considering the composition dependence of the 
free energy density. 

To find the divergent part of Eq. (|189[) for Cij (k) , we 
will need high-g asymptotic expansions for the screened 
interaction Gy(q+), and the function Di(k, q , — q+). 
To the required order in an expansion in powers of 1/q, 



Gy(q±) 


= G<°> -+ 


G (o) 


= V --f 
Y1V 


G<±> 


12; 


G« 


- - ( 

2Z 2 V 




- 2U Z 
- 



(±)+gW+g 



(x) 



iVi 



2„,2 



few 
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To obtain a corresponding asymptotic expansion of Di in 
the high-g limit, for fixed k, we note that 



A(k,q) 



(3) 



(K, Q , -Q+) 



9(Kf 



(194) 



and use expansion (| 1T2[) of uj( 3 \ This yields 



A(k,q) =s £>. 



(0) 



(i) 



D 



(o) _ 



(i) 



12Zj 
v<7 2 
3k 



(195) 



/c 2 (3 + 4a 2 ) 



72L 2 



v 2 N i9 (Kf) 



where AT 2 = k 2 Nib 2 /6 = k 2 N z v/(6k). The key simplify- 
ing feature of this expansion is the fact that the leading 
order contribution to Z5j(k, q) is independent of k. 

By substituting these expressions into Eq. (| 1 89[) . and 
keeping only terms of 0(1) and 0(l/q 2 ) in the integrand 
that lead to UV divergent integrals, we obtain a UV di- 
vergent contribution to SCij (k) as a sum: 



where 



sc$> 



self 



\ f D^D^C^C^ 



D (0) D (0) G (x) G (0) 



SC. 



(i) 



(196) 



(197) 



SC. 



D W D f)[G(+)&-)+2&VGM} 



G (0) G (0) 



Evaluating the integrals with respect to q yields 



SC 



(0) 



sc 



(x) 

ij 

.(1) 



lih v a 3 



12tt 2 Z 2 



12i> LL 



<5G}; ; = 



7k%ljV 2 , 
12tt 2 P 



(198) 



3liljV 



7T 2 Z 3 V #1 



No 



+ 



A 



2tt 2 Z 2 liV 4ff (^ 2 ) A^(K 2 ) 



A 



The divergent part of the corresponding apparent x pa- 
rameter is a sum of the three terms 



*x (0) 



(h~h) 2 v A 3 
8n 2 l 2 3 
6Z 2 Z| 

XoA 



7T 2 / 3 

7fc 2 (Z x 



Z 2 ) 2 ^ 



24tt 2 ; 2 
3(ii - / 2 ) 2 

2tt 2 Z 3 
3(Zi - Z 2 ) 

2tt 2 [ 2 



A 

'2 



(199) 



A, 



A 



/ 2 
'i 



A 25 (^ 2 ) N l9 {K{) 



A 



Note that both \^ an d X vanish in the case l\ = I2 
of equal statistical segment lengths. 

The UV convergent contribution <5G Z * (k) is given by 
the difference between the exact integral expression of 
Eq. (|189|) and the sum of the UV divergent inte- 
grals given in Eqs. (|196|) and (|197|) . The resulting 
UV convergent integral can be non-dimensionalized by 



25 



arguments similar to those used to non-dimensionalize 
5il*Ak), which in this case yield a non-dimensional inte- 
gral times a prefactor of v / (N N 1 / 2 ). The corresponding 
non-dimensional form of S\*(k) is given in Eq. (|55[) . 



XI. CORRELATIONS IN DIBLOCK 
COPOLYMER MELTS 

We now consider the calculation of SSij (q) for a diblock 
copolymer melt. The calculation is closely analogous to 
that given above for a blend. 



A. Two-Point Functions 

The two-point intramolecular function f2jj(q) for a 
Gaussian diblock copolymer of length N with blocks of 
length fiN and /2-W, is a matrix 



N 

^y'(q) = — ffy(q) 



(200) 



subsection IV Bl for a binary homopolymer blend. The 
UV divergent contribution is of the form 



Sf ~ Sf (0} + Sf (x> + Sf 



ix) 



c(l) 



(203) 



The expressions for <5/(°) and Sf^ are identical to those 
given in Eqs. ([52")) and ([M]) . respectively, for a blend 
of the same composition. As for a homopolymer blend, 
the sum of these two terms yields the quantity 5f\ oca i — 
Sf(o) + 8 fix). The remaining 0(1/ N) contribution is 
given in a diblock melt by 



Sf 



(i) 



3 l\ — I1I2 + I2 ^ 



2ir 2 vl 



(i) 



Here, 



' (2) 

contribution to fi:/ 



Ea^iVl where 0%'"' is the 0(N/Q 4 ) 



N 



(2,1) 



(204) 



(q). 



where 



ffi eie 2 
eie 2 32 

2( e -A<3? 

-/< 



(1 



1 + fiQD/Qt 
)/Ql 
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with Q 2 = q 2 Nb 2 /6. The high-q asymptotic behavior 
of these function may be obtained by dropping all terms 
that contain factors of e - -'^ . 

The high-g behavior of the propagator Gy (q) may be 
approximated to the required accuracy by an expansion 

[ll-hh + lpv , 2v 

+ Wn + P x " z%3 

Here, I and Z^ have the same values as in a homopolymer 
blend of the same composition, with <pi — fi. That is, 



I 



Z = 



flh + /2^2 

_ /2^2 hhhh 

hhhh 



-mi 



(202) 



The only difference between this expansion and the cor- 
responding expansion of Gy (q) in a homopolymer blend 

of equal composition is the nature of the term Gy (k) 
that is proportional to 1/N and independent of x, which 
is the first term in the second line of Eq. (|202p . 



B. Free Energy Density 

The one-loop contribution Sf to the free energy density 
of a homogeneous diblock copolymer melt can be calcu- 
lated by a procedure closely analogous to that given in 



Intramolecular Correlations 



In the diblock case, we calculate 8uj\ ■ (k) using Eq. 
(|15ip . To analyze the divergence of this expression 
we need to identify the high-q behavior of its compo- 
nents. As in the blend case, this integral has a diver- 
gent part proportional to A, which may be obtained by 
using the dominant contributions to r/>yjy(k, — k, q, — q) 

~ ( 2) 

andG^(k). As in a blend, the dominant contribution to 

G kl (k) is independent of the values of the indices, and 
is given by Eq. 

The high-g behavior of V'yjw ^ s analyzed in the ap- 
pendix. As for the case of the blend, we find a 
simple expression for the dominant asymptotic behav- 
ior ip^ which largely determines the form of the 
result: The dominant contribution to the function 
^jy(k, — k,q, — q) arises only for elements with k = I, 
and is given by 



-k,q,-q) 



288a 2 ft5g>(k) 
d(b 2 ) 



q% 2 



(205) 



For i = j, the derivative with respect to b 2 . is nonzero only 
for i = j = k, while for i =/= j, the derivative is nonzero 
for both k = 1 and k = 2. The dominant contributions to 
elements of tp^jli (k, — k, q, — q) with k ^ I are all found 

to be 0(1/ q 6 ) or smaller, and thus do not contribute to 

f 21 

the divergent part of Eq. (|15ip for SQ^ (k) . 

Using this asymptotic result for ip^ 4 \ the UV divergent 
contribution to Eq. (|15ip can be written as a sum 



M, 2) (k) 



E 

k 



d(b[ 



(206) 
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where S(b 2 ) is again given by Eq. (| 1 8T[) . For the diblock 
melt, as for the blend, we thus find that the divergent 

(2) 

contributions to Slj^ (k) can be absorbed into a renor- 
malization of statistical segment lengths. 



where k ^ j. In order to calculate the one-loop contri- 
bution to 5Cij(h), we will need a quantity 



A (1; (k,q) 



^ ^ikl 
kl 



(k,q) 



(214) 



D. Direct Correlation Function 



This is given by 



The one-loop contribution to the direct correlation 
function in a diblock copolymer melt can be written as 



JCy(k) 



1 



D ikm (k, q)G fci (q + )Dfil (k, q)G mn (q_ 



where, in this context, we define 



(207) 



A^k,q) = f\^g(k,q_,-q + ) , (208) 
with Diki(k, q) = Diki{— k, — q). As for the binary 



S(3) 



blend, we expand the three point function Qj fci — 

(209) 



fl^(k, q_, — q_)_) in the high-g limit as a sum 



m = & 3 k f + 

jkl jkl 1 jkl 

and Diki (k, q) as a corresponding sum 

Aw ~ L» 4fc; + L> ifcJ 

where f§ 0) and d|°J are 0(l/q 2 ) and Clf k f and D$ 
are C(l/g 4 ). 

A detailed analysis of the high-g behavior of 0,j kl , 
which is outlined in the appendix, shows that the leading 

order contribution fi^y , which is of 0(1/ q 2 ), is nonzero 

only for k — /, and has four elements &^ k that are all 
of the form 



(210) 



^(k,q_,-q + ) = Q^(k)^ (211! 



- f)( 2 )/ 



121* 



It follows from definition Eq. (|208[) and this expression 
that the corresponding leading order contribution to D^i 
is 



D^ l (k,d)=6 ik 6 kl D^(c l ) 



(0), 



where 



A (0) (q)-^ 

vq z 



(212) 



(213) 



is independent of k, and is the same function as that 
found in the binary blend. 

The only nonzero elements of the subdominant contri- 
i(3,i) _ o(3,l). 



bution Cl\ k \ = Q\ k \ '(k, q_,-q+) are: 
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^ 2 (3 + 4a 2 ) 5 ,-4/,] 
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3 -e,\K 2 (3 + W)e k -^ 
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36ljlkej 1 
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= 



v fc2(3+4a2) -^^ 



71 ( 215 ) 



where we have defined a vector Li with components 



Li = fill + eill - eihl 2 
L 2 = hi 2 



2 T S2I1 ~ 62/1/2 



(216) 



The expression for in a homopolymer blend given in 
Eq. (|85[) may be recovered by the replacements Li — > /? 
and g^ 1 /N -> %/( 5 iiVi). 

When expressed in terms of and the quantity -Di 
defined above, the integral expression for the divergent 
part of <5Cjj(k) in a diblock copolymer melt is identi- 
cal to that given in Eq. (|197j) for a homopolymer blend. 
The only differences between the expressions obtained for 
5Cij in a diblock copolymer melt and that in a homopoly- 
mer blend of the same composition arise from the use of 
different expressions for the the fc-independent part of 
G?y [given for a diblock by the first term of the second 

line of Eq. [202], and for Z)f } [given by Eq. ([215]) ]. More- 
over, the integral expressions for 5C^' (k) and <5Cy (k) 
are identical to those obtained for a homopolymer blend 
of the same composition. Only the contribution SCj-j (k) 
differs from that obtained for a corresponding blend. 
This is given by 



SC 



(i) 



7k 2 kl 3 v 2 

12tt 2 /" 2 
3vklj(ll 



A 



hh 



3± A 



3v 



2Ntt 2 V 



TT 2 Nl 3 

-Xg^Lklj+hLkg^hk (217) 



The first line, the term proportional to fc 2 , is identical 
to the corresponding expression for a blend of the same 
composition. The second line, which arises from the inte- 
gral involving the fc-independent part of G", is different 
because of the use of different expression for & x >. The 
third line is different because of the use of a different 
expression for . 

The corresponding expression for the UV divergent 
contribution to <5x(k) is similar to that obtained for a 
blend of the same composition. The expressions for Sx^ 
and Sx^ x ' terms are identical to those obtained in a ho- 
mopolymer blend. The expression for Sx^ (k) is differ- 
ent, but is retains the property that 8x^ = in the limit 
Zi = h of equal statistical segment lengths. 
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XII. END EFFECTS 



where 



We show in this section that the form of our results for 
the UV divergent part of <5CV, (k) are consistent with the 
existence of a UV divergent one-loop contribution to the 
free energy of the form proposed in Eq. (|56|) . Here, it is 
convenient to start from the more general expression 

SF int = J dr[ S /local + ^2 da^a + -DijVCiVCj } , 

(218) 

in which the gradient-squared term is expressed in terms 
of monomer concentrations, as would be required to de- 
scribe a slightly compressible liquid. The postulated free 
energy is the sum of a Ginzburg-Landau like functional 
of the monomer concentrations plus additional free ener- 
gies arising from chain ends and (for a diblock) from the 
junction that connects the two blocks. 

The excess free energies arising from chain ends and 
junctions are assumed to depend linearly on the concen- 
trations of these defects because the defect concentra- 
tions are dilute. Direct defect-defect interactions, which 
would yield contributions to S^ 1 (k) of 0(l/N 2 ), could 
become important in systems of relatively short chains 
with strongly interacting end-groups, but do not appear 
in our model within the one-loop approximation. 



A. Monomer and Defect Chemical Potentials 

We compare our one-loop results to a form of SCFT 
in which the average monomer concentrations are cal- 
culated from those of ideal gas reference system with a 
self-consistent field Hamiltonian 



U = t/chain ~ >^ hj * Cj + ^ d a * 1p a 



(219) 



in which a chain end of type a = 1 or a = 2, or a junction 
(a = J) at position r is penalized by a free energy ip a (r). 
In the absence of any external field, monomers of type i 
are subjected to a field 



hi = Uij * cj + 



5(5^ 



intJ 



5c, 



(220) 



in which 5F- ln t is given by Eq. (|218p , and S /Sc-i represents 
a functional derivative. 

To construct an RPA calculation of S(k), we consider 
the deviations from a homogeneous reference state in- 
duced by a small external field hi that couples only to 
the monomer density. This external perturbation will 
induce deviations 5ci and 5d a in monomer and defect 
concentrations, respectively, and deviations Shi an d Sip a 
in the conjugate fields. These deviations in chemical po- 
tential fields are given, to linear order in the deviations 
in the concentrations, by 



5hi = U-j * Scj 



+ Vi a * Sd a - hi 



Vu 



6c,. - UlJ 



5 2 (6F mt ) 
SciScj 



frtpa 

Sd 



(222) 



Here, U[j is an effective monomer-monomer interaction, 
and Vi a where is an effective monomer-defect interaction. 
Here, and throughout this section, we use '*' to represent 
spatial convolution only, and display the monomer and 
defect type indices explicitly. 

In Fourier space, the effective interaction ?7/-(k) is 
given by 



E7f,.(k) = ^(k) 
where 



d 2 {6 /local) 

dcidcj 



dcidcj 



(223) 



(224) 



Here d a is the defect density in the homogeneous ref- 
erence state: For hopolymers, the ends densities are 
d a = 2cf>a/(vN a ); for diblock copolymers, both ends and 
junction densities are d a — l/(vN). All derivatives are 
evaluated in this homogeneous reference state. Note that 
our definition of UL as a second functional derivative of 
the interaction free energy yields a contribution Dijk 2 
that arises from the (postulated) gradient-squared con- 
tribution to (5 -Pint ■ 



B. Generalized RPA 

To complete the RPA linear response calculation, we 
must combine the above with a description of the linear 
response of an ideal gas of polymers. The linear response 
of the monomer and defect concentrations to deviations 
in the above combination of fields is given by 



Sci 
Sd a 



5hj 



R aj * 5hj 



Rip *Sip 
- E a * Sipp 



(225) 



Here, f2y is the ideal-gas intramolecular correlation be- 
tween monomers of types i and j, Rj a is the intramolec- 
ular correlation between a defects and j monomers, and 
E a p is the intramolecular correlation between defects of 
type a and (3. All of these functions are diagonal in a ho- 
mopolymer blend (i.e., are nonzero only for i = j , i = a 
or a = /3), and all become non-diagonal in a diblock 
copolymer melt. 

By combining linear response equations (|225[) with self- 
consistency conditions (|221[) , it is straightforward to show 
that, in a Fourier representation, #Cj(k) is given by 



Ho 



(221) 



Sc t = Ai 3 [-B jk 5c k + hj] 



(226) 
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where all quantities are implicitly functions of a 
wavenumber k, and 



Bjk = Uj k - Vj a E a pVj. 



iik 



(227) 



Solving for the inverse response/correlation function 
S^ik) = 6h j (k)/8c i (k) yields 



s ij - Aj 1 
= ft 1 



Ui,. 



in which 



6SJ = 



d{8h 



local J 



deidej 



D l3 k 2 



SC. 



(228) 



(229) 



is the correction to mean-field theory, and where SC^ 
is a contribution arising from the interaction of the end 
and junction 'defects' with the monomer concentration. 
This quantity is given by a sum 



-sd? 1} - 5C {d2) 



-SC, 



(d2) 



Vic 



VaPak^ki + ^ik R kaV aj 
T 



P^k^H R W - E > 



Note that the quantity SC^ 1 ' depends linearly upon the 

effective interaction Vi a , while SC^ 2 ^ is second-order in 
V ia . 



a/3 



V, 



(230) 



C. Comparison to One-Loop Results 



In both blends and diblock copolymer melts, the re- 
sult from the one loop approximation is thus identical 
to that obtained from the generalized RPA, except for 
the absence in our one-loop results of the contribution 
SC[ d predicted by RPA result. Note that the contri- 
bution SC^j that is "missing" from the results of the 
one-loop approximation has a qualitatively different de- 
pendence on both Vi a and A than terms that are found 
in the one-loop calculation: The quantity 5C[ 2d ^ is sec- 
ond order in the strength of the coupling Vi a , and so 
would yield a contribution SC\ 2d ^ oc A 2 in a theory in 
which Vi a oc A, while the terms that do appear in our 
one-loop results are linear in Via, and proportional to A. 
The one-loop calculation of SCij yields only terms pro- 
portional to A and A 3 , so it appears that the missing 
contribution predicted by the RPA could not possibly be 
generated by a one-loop approximation. We assume that 
the discrepancy occurs because terms that are quadratic 
in the magnitude of the one- loop contributions to Vi a ap- 
pear only at second order in a loop expansion of S (k). 
The A-dependence of the missing contribution appears to 
be consistent with this conjecture: The power counting 
analysis of the loop expansion given in Sec. IXIIII and in 
Ref. [13] indicates that the dominant UV divergence of 
the two-loop contribution to SC will be C(A 4 ), to which 
we expect to find subdominant corrections of 0(A 2 ) and 
O(lnA). While this conjecture could be proved only by 
analyzing the two-loop theory, it seems clear that the 
missing terms cannot be produced by a one-loop approx- 
imation, and so may be neglected when making this com- 
parison. 

Below, we present some details of the application of the 
generalized RPA to binary blends and diblock copolymer 
melts, respectively. 



To test our phenomenological model, we compare our 
explicit one loop calculations for the SCffl (k) to the re- 
sults of the above generalized RPA calculation, while us- 
ing the explicit expressions for tp a given in Eqs. (|59p and 
(l61j) . and treating Dij as a free parameter. We find that, 
in both blends and diblock copolymer melts, the quantity 
SC. 



(i) 



( obtained in the one-loop approximation for SCij (k) 
can be written as a sum of the form 



D l3 k z 



SC { x dl \k) 
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where SC\ d1 ^ is the first-order defect contribution given 



in Eq. pgDj) . and 



Di 



V 2 klj 



2/2 



A 



(232) 



The same value is obtained for in blends and diblock 
copolymer melts of the same composition. The scalar 
coefficient given in Eq. (|62p is obtained by requiring that 
Vci = — Vc2 in an incompressible liquid. 



1. Binary Blends 



For a binary blend of homopolymers, 
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where gi = g{K 2 ) is the Debye function, and where ti 
(1 — e~ Ai )/ K 2 for a homopolymer. This yields 



SC. 



(di) 



27T 2 /" 2 



• v ,.'/. ^ j a, 



A (234) 
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By using the identity = 1 — g^K 2 /2 and the definition 
of K 2 = k 2 v/(6li), we obtain the alternative expression 



6C™ = Wij 



2n 2 l 2 



u 



By substituting this expression for SC[j into Eq. (|23ip . 
and using Eq. (|232|) for Dij, we reproduce the one- loop 
result for 5C\f given in Eq. (TTggj) , 



A 



47r 2 r 2 
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3. Diblock Copolymer Melts 

For a diblock copolymer melt, we obtain 
3uL 



V 7 , 
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(ii-^rA 
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SijV 1 Ng i] 
v~ 1 e i (a = i) 



V Ci 



(a ^ i, J) 



(237) 



where gij, gi, ej are defined for a diblock copolymer by 
Eq. PJT]| . This yields 



5Clf = Wij 



3vLl 



2w 2 Nl 2 



^-etlj + ke k g~l\ A (238) 



By using the identity ej = /j — giKf/2, we obtain the 
alternative expression 



A-^A 



9 l k L kh + hL k g k j j- 4 ^ 2p 

(239) 

where Lfe is given by Eq. (|2 1 6[) . By using this expres- 



sion for <5Q" ' and Eq. J232|) for D l3 in Eq. (|23l]) . we 
reproduce Eq. (f2T7) . 



XIII. BEYOND ONE-LOOP 

In this section, we briefly look beyond the one-loop the- 
ory. We consider the structure of UV divergences of an 
arbitrary diagram in the unrenormalized diagrammatic 
perturbation theory, at the level of naive power counting. 
We also consider the structure of a renormalized loop ex- 
pansion, and argue that this should yield an asymptotic 
expansion of corrections to SCFT in powers of 1/VN. 

For this purpose, we use the diagrammatic formalism 
presented in Ref. [13]. There, it was shown that cor- 
rections to the Gaussian or one-loop approximation for 
the grand potential In Z for a polymer liquid could be 



expressed as an infinite sum of connected diagrams of 
— G bonds and SI vertices, with no root sites. Here, we 
consider the dominant UV divergence of an arbitrary di- 
agram in the expansion of the free energy density In Z/V 
for a model with a wavenumber cutoff A. For generality, 
we consider a model of continuous Gaussian polymers in 
a space of arbitrary spatial dimension d. We consider a 
generic connected diagram with B bonds, each represent- 
ing a factor of — G, connecting V vertices, each represent- 
ing a factor of Q( n \ where n is the number of attached 
bonds. The number L of loops, or independent wavevec- 
tor integrals, in such a diagram is L = V — B + 1. For 
simplicity, we consider only the dominant UV divergent 
to the free energy density of a one-component liquid (or 
of a blend with \ = and b\ = 62), without attempting 
to analyze the dependence of the free energy density on 
X, or on the composition of a blend. 

The analyze the UV divergence of an arbitrary dia- 
gram, we must characterize the asymptotic high-g be- 
havior of G(q), and of Sl' n )(qi, . . . ,q„) for arbitrary n. 
The dominant g-dependence of the screened interaction 
G(q), given in Eq. fl57D, is 



G(q) ~ {qbf/c 



(240) 



where c = 1/v is the monomer concentration. The func- 
tion can be expressed, for a Gaussian homopoly- 
mer in any spatial dimension d, as a product of the form 
n(")( qi ,...,q„) = cN n - 1 u^{Q l ,...,Qn), where c/N 
is the number concentration of polymers, Qi = qi-R, 
and R = by/ N/6. To characterize the high-g limit, it 
is useful to consider the limit of infinitely long chains, 
N — > 00, at constant monomer concentration c. In 
this limit, all nonzero wavevectors become "large" com- 
pared to 1/R. The function f2(")(qi, . . . ,q n ) generally 
approaches a nonzero limit as N — > 00. For QS n > to ap- 
proach a value that is independent of N in this limit, 
i!)W(Qi, . . . , Q n R) must approach a homogeneous scal- 
ing function 

lim a;W(AQ 1) ...,AQ„)=A 2 ( n - 1 ^W(Qi,...,Q n ) 

iV— 'CO 

(241) 

when |Q| 3> 1 for all n arguments. At a power count- 
ing level, this implies that the function u)(") is of order 
q 2 (™- 1 ) w hen all of its arguments are of order Q. Corre- 
sponding, Qy 1 ' must be of order 



(qb) 2n 



(242) 



when all of its arguments are of order q ^> 1/R. 

To count powers of q in the Fourier integral associated 
with an arbitrary diagram in the expansion of In Z/V , 
it useful to associated one of the n factors of (qb) 2 in 
the denominator of Eq. (|242[) for &S n ' with one end of 
one of the n bonds that must be attached to the asso- 
ciated vertex. This method of counting counting leaves 
an overall factor of c(qb) 2 for each vertex, and a factor 
of c~ 1 (qb)^ 2 for each bond in an incompressible liquid. 
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In a diagram with L loops, in d dimensions, we obtain 
an additional factor of order q dL from the integration 
over L wavevectors. By combining these factors, we find 
that the contribution to In Z/V from a generic connected 
L-loop diagram with a cutoff wavenumber A is of order 

A d (A d - 2 c/6 2 ) i - 1 (243) 

In the physically relevant case d — 3, we may define a 
packing length p = l/(cb 2 ), to rewrite this as 

A 3 {A P ) L - 1 (244) 

The order of UV divergence is thus expected to increases 
by one factor of Ap at each order in the loop expansion. 

In the one-loop theory considered here, with d = 3 and 
L = 1, this argument yields a UV divergent contribution 
of order A 3 . This agrees with the results of our explicit 
one-loop calculation of free energy density of an incom- 
pressible liquid. Corresponding L-loop contributions to 
the functional derivatives of F [(c)], such as the second 
derivative S' _1 (fc), are expected to have the same depen- 
dence on A as the underlying free energy contributions, at 
eac order in the loop expansion. Thus, for example, the 
dominant UV divergent one-loop contribution to S^ 1 (k) 
in our explicit one- loop calculation is also of order A 3 . 

We have shown that this UV divergence of the one- 
loop theory can be removed by renormalization. Let 
us assume, for the moment, that this procedure can be 
extended to arbitrary order in a loop expansion. We 
imagine that the contribution of an L-loop diagram to 
the renormalized perturbation theory can be obtained 
by subtracting an asymptotic approximation for the inte- 
grand of the corresponding Fourier integral that is accu- 
rate at high wavenumbers, q 3> 1/R. This will generally 
leave a UV convergent contribution that arises primar- 
ily from wavenumbers of order 1/R, due to deviations of 
the integrand from the asymptotic approximation at low 
wavenumbers. The resulting contribution to the renor- 
malized perturbation theory is expected to be similar in 
magnitude to the value of the unrenormalized integral 
evaluated with a cutoff A ~ 1/R. This suggests that 
the renormalized perturbation theory (if one exists) will 
yield a loop expansion in which the L-loop correction to 
the SCFT free energy functional will have a prefactor of 
order 

i? 3 (1) ~ vNN L / 2 (245) 

We have used the fact that R/p = N 1 / 2 b z /v = N 1 / 2 to 
obtain the second expression in the above. The value 
of each diagram will be given by this prefactor times a 
non-dimensionalized convergent integral whose value is 
a dimensionless function of %N and the other SCFT di- 
mensionless paramters. If our renormalization procedure 
can be extended beyond the one-loop level, the resulting 
renormalized loop expansion is thus expected to yield 
an asymptotic expansion of the free energy in powers of 

i/Vn. 



Eq. (|245p may also be obtained by dimensional anal- 
ysis: If we non-dimensionalize all lengths in units of 
R = ^/Nb in the functional Taylor expansion of the sta- 
tistical weight L that appears in the Edwards functional 
integral, we obtain an expression for L as the product of 
a non-dimensionalized functional times a large parameter 
V~N. The existence of such a large prefactor is known to 
imply 2 i that the loop expansion (or, more precisely, the 
renormalized loop expansion) will yield an expansion in 
powers of N~ x / 2 . 

XIV. POLYMER SOLUTIONS 

UV divergences also appear in the excluded volume 
problem for polymers in good solvent. The diagram- 
matic formalism and power counting arguments given 
above for an incompressible liquid can also be applied 
to the excluded volume problem. A comparison of the 
two problems in the same language is instructive. 

In the standard Edwards mode l 11 ' 12 of a polymer in 
good solvent, the two-body interaction is approximated 
by a point-like effective interaction U[r) = aS(r), or 
U(q) = a, where a is an effective excluded volume. The 
one-loop theory given here for incompressible liquids is 
closely analogous to Edwards' one- loop theory of solu- 
tions. The screened interaction in the original Edwards 
theory is given, in the limit N — » oo, by 

G ~ 1 (?) = ^|+«~ 1 ( 24 6) 

where the first term is the limit qR 3> 1 of fl (q). In this 
theory, we thus obtain an essentially unscreened inter- 
action G(q) ~ a for fc£ 3> 1, and a screened interaction 
identical to that given in Eq. (|240|) for fc£ -C 1, where 
£ ~ b/y/ca is the Edwards screening length. 

The one-loop theory for a nominally incompressible liq- 
uid is thus very similar to the Edwards theory with a 
cutoff wavenumber A < In studies of non-dilute 

polymer solutions one is generally interested in the de- 
pendence of the radius of gyration, osmotic pressure etc. 
upon contributions from wavelengths less than £, which 
determine the concentration dependence of these quan- 
tities. To retain this information, one must use a cutoff 

A»r 1 - 

The UV divergence contributions from the regime 
q£ 3> I may be analyzed by repeating the power counting 
arguments given above for an incompressible liquid, while 
using the unscreened interaction G(q) = a, rather than 
the screened interaction G(q) cx q 2 . By this method, we 
find that the contribution to the free energy density of 
an L-loop diagram with V vertices is of order 

A d (A6 )(i-i)(d-4)-2y (a/6( i ) L-i (ac) y (247) 

In the infinite dilution limit, in which one considers only 
the interaction among monomers of a single chain, the 
only relevant diagrams have only one vertex, V = 1. In 
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this limit, the above expression becomes equivalent to the 
known result 28 for the naive degree of divergence of dia- 
grams in the perturbation theory for swelling of a single 
self- avoiding chain. 

In the physically relevant case d = 3, the only UV 
divergent diagrams in the expansion of the free energy 
density for a solution are the one-loop diagram with one- 
vertex (L = V = 1), which diverges as A, and the two- 
loop contribution with one vertex (L = 2 and V = 1), 
which exhibits a logarithmic divergence in d = 3. The 
strongly divergent one-loop diagram was correctly iden- 
tified by Edwards with a divergence in the free energy 
per monomer due to interactions between pairs of nearby 
monomers on the same chain. Edwards removed this di- 
vergence by subtracting the free energy per monomer of 
a single isolated chain from the total free energy. All of 
the diagrams involving more than one vertex, which are 
necessary to calculate, e.g., the second virial coefficient, 
are UV convergent in d = 3. 

Alexander-Katz et al2& have recently considered the 
UV divergence of the chemical potential in a stochastic 
field-theoretic simulation of a polymer solutions. In this 
context, the UV divergence shows up as a dependence 
of the polymer chemical potential \x upon a spatial dis- 
cretization length Ax. They found that that the UV 
divergence of /x found in their simulation could be re- 
moved by subtracting a one-loop approximation for the 
free energy per monomer of a system of non-interacting 
chains from their simulation results. (The above analysis 
suggests that a remaining logarithmic divergence should 
have remained, but they reported no evidence that sug- 
gests this). Our analysis of the incompressible liquid 
model indicates that this simple subtraction will not be 
sufficient to remove the UV divergences from analogous 
simulations of dense polymer mixtures, in which the di- 
vergence appears at all orders in the loop expansion. 

The appearance of negative exponents of A in all but 
a few diagrams in the expansion of the dilute solution 
free energy in d = 3 is a symptom of the fact that most 
diagrams in this theory are infrared (IR) rather than UV 
divergent. In the infinite dilution limit, in which only 
diagrams with V = 1 are relevant, the IR divergence is 
cutoff at a wavelength of order R ~ yNb. In this limit, a 
generic diagram diverges with increasing chain length N 
as N d / 2 N { ^ d ^ L ^/ 2 for all d less than the upper critical 
dimension d c = 4. In a semidilute or concentrated solu- 
tion, this IR divergence is cutoff at the screening length 
£. The quantities that we calculate in this paper, other 
than the free energy density (i.e., the direct correlation 
function, which is analogous to the second virial coeffi- 
cient, and the interaction-induced change in the single- 
chain correlation function) are actually UV convergent 
and IR divergent in the dilute solution problem. It is, of 
course, the IR divergence of the excluded volume prob- 
lem that makes the problem interesting, and that leads 
to nontrivial scaling behavior for a self-avoiding walk. 

The above analysis makes it clear that the excluded 
volume problem is "renormalizable" in the sense that this 



word is normally used in quantum or statistical field the- 
ory, and that the theory of incompressible liquids consid- 
ered here is not. A field theory is normally said to have 
a renormalizable UV divergence if the divergence can be 
absorbed into a finite number of measurable parameters, 
such as the mass and charge in quantum electrodynam- 
ics. More precisely, it is usually required that only a 
finite number of vertex functions contain a primitive UV 
divergence in spatial dimensions d less than or equal to 
an upper critical dimension d c £^ For this to be true, it 
must generally be the case that the degree of UV di- 
vergence of all vertex functions decrease with increasing 
order in a loop expansion for all d < d c . For any d < d c 
only a few low order diagrams thus remain UV diver- 
gent, but the theory instead becomes susceptible to the 
appearance of IR divergences in (in different contexts) 
the limit T — > T c or N — > oo, leading to nontrivial crit- 
ical phenomena. The Edwards model for polymers in 
good solvent fits this description, with an upper critical 
dimension d c — 4. In the model of incompressible liquids 
considered here, however, we see from Eq. (|244j) that the 
degree of UV divergence increases with increasing order 
in a loop expansion for all d > 2, implying that d c = 2. 
This model would thus normally be said be unrenormal- 
izable in d = 3. The difference between the upper critical 
dimensions of the two theories is a result of the replace- 
ment of the unscreened point-like interaction, which is 
independent of wavenumber, by a strongly wavenumbcr 
dependent screened interaction, G(q) cx q 2 . 

The notion of "renormalizability" used in this paper is 
thus different from its usual meaning in field theory, and 
is specific to the physics of dense polymer mixtures. In ei- 
ther context, a theory is said to be renormalizable if all of 
its UV divergences can be absorbed into the parameters 
of an appropriate phenomenological model. In addition, 
for renormalization to be useful, it is required in both 
contexts that the functional form of a phenomenological 
model that is sufficiently flexible to absorb all UV diver- 
gences also be sufficiently constrained to allow nontrivial 
predictions to be made. A theory is properly described as 
unrenormalizable if the functional form of the theory re- 
quired to absorb all UV divergent parts of the free energy 
functional is so flexible that it could describe all conceiv- 
able behavior. In the present context, the appropriate 
phenomenological model is a SCFT of Gaussian chains. 
The only constraint that we place on the functional form 
of this theory is that interaction free energy functional 
(excluding the contributions of chain ends and junctions) 
be independent of chain length N and chain architecture. 
We show here that, to first order in a loop expansion, a 
SCFT of this form is recovered as the N — > oo limit of the 
true free energy functional, and that the sensitivity of the 
theory to the cutoff (or more generally, to monomer scale 
chemical details) can be absorbed into this N — > oo limit. 
The renormalized perturbation theory makes nontrivial 
predictions about iV-dependent, UV convergent correc- 
tions to this theory. If the same ideas were applied to a 
liquid of point particles, there would be no constraint on 
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the functional form of the free energy, and so the theory 
would have no predictive power. It is thus the appear- 
ance of 1/N as a small parameter that makes it possible 
to extract useful information from a theory that would 
be considered unrenormalizable by the usual rules of field 
theory. 



XV. CONCLUSIONS 

Several previous calculations have shown that predic- 
tions of a coarse-grained model of polymer liquids, when 
extended beyond the mean-field level, depend very sensi- 
tively on the value chosen for an arbitrary coarse-graining 
length. In light of this, it is reasonable to ask whether 
such models can make any unambiguous physical pre- 
dictions. We know how to extract physical predictions 
from a UV divergent theory only if we can absorb all 
UV divergences into the values of a finite number of 
phenomenological coefficients. The question of whether 
such coarse-grained theories can predict anything thus in- 
evitably boils down to the question of whether they are, 
in some sense, renormalizable. The definition of "renor- 
malizability" that we introduce here (as discussed above) 
is based upon the assumption that a renormalized SCFT 
with an interaction free energy that is independent of N 
is obtained as the limit N — > oo of the true free energy, 
and that all UV divergences can be absorbed into this 
asymptotic theory. 

In this paper, we have explicitly demonstrated the 
renormalizability of the one-loop approximation for 
S (k) in the auxiliary field approach, for arbitrary k, for 
both polymer blends and diblock copolymer melts. We 
showed that all UV divergent contributions to S^ 1 (k) 
can be absorbed into changes in the values of the pa- 
rameters of an appropriate form of SCFT. To make this 
work, it was necessary to allow for renormalization of all 
of the parameters of the standard SCFT, i.e., of the sta- 
tistical segment lengths as well as the local interaction 
free energy. We also found that it was necessary to al- 
low for some free energy contributions that are plausible 
on physical grounds, but not usually considered, i.e., a 
square-gradient interaction free energy and excess ener- 
gies for chain end and for junctions in block copolymers. 

At a conceptual level, this analysis is important be- 
cause it provides evidence for the logical consistency 
of the study of coarse-grained models, and for our as- 
sumption that some form of SCFT is exact in the limit 
N — > oo. Because we have analyzed only one vertex func- 
tion, S^ 1 (k), only to first order in a loop expansion, our 
calculation provides a consistency check, but not a proof 
of renormalizability. This is quite different from the sit- 
uation in quantum and statistical field theory, in which 
methods were developed long ago to prove renormaliz- 
ability to all orders in perturbation theory. An analo- 
gous proof would require an enormous generalization of 
the explicit analysis given here. 

At a practical level, our analysis is important as a nec- 



essary step in the development of a rigorous renormalized 
perturbation theory of corrections to SCFT: Identifica- 
tion and removal of all UV divergences is a prerequisite 
to the systematic study of the long wavelength physics 
that coarse-grained models are intended to describe. In 
subsequent work on this subject, we will focus on exam- 
ination of physical predictions of the renormalized one- 
loop theory for corrections to SCFT, and on extending 
the theory beyond the one-loop level. 

Our procedure for extracting the UV convergent con- 
tributions of physical interest is, for the moment, to sim- 
ply subtract the results of our analytic calculation of the 
UV divergent contribution to each quantity from the re- 
sults of a numerical evaluation of the underlying Fourier 
integral, using the same finite cutoff wavenumber A in 
both calculations. To make this procedure work, we must 
identify and subtract all UV divergent contributions to 
the unrenormalizcd integrals, including the 0(A/N) di- 
vergences discussed in Sec. IXIII We have confirmed 
numerically that this procedure yields results that are 
nearly independent of A for large values of A, and that 
converge in the limit A — > oo. It is worth noting that this 
numerical procedure provides a very stringent test of the 
correctness of our results: Any error in cither our ana- 
lytic calculation of UV divergent contributions or in our 
numerical integration would destroy the required cancel- 
lation of UV divergent terms. 



APPENDIX A: SINGLE-CHAIN CORRELATIONS 

In this appendix, we provide further details of the 
derivation of asymptotic high-g expansions of the three- 
and four-point correlation functions for both homopoly- 
mers and diblock copolymers. The leading order terms in 
the required expansions were obtained for homopolymers 
in Sec. |IX] 



1. Generalization to Block Copolymers 

To begin, we generalize to block copolymers the discus- 
sion of multi-point correlations for a Gaussian homopoly- 
mer that was given in subsection IIX Al We consider a 
n-point correlation function 

JV 

^,i 3 (ki,...,k n ) = |cr S ( e ^' R fe)) , (Al) 
o 

Here, d n s should be understood to indicate an integral 
over all values of s\, . . . , s„, subject only to the constraint 
that, in a block copolymer, the integral over Sj must be 
taken over only the block of the polymer that contains 
monomers of type ij. (For simplicity, we will only dis- 
cuss situations in which each block is chemically distinct). 
For example, in a diblock copolymer for which monomers 
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with < s < fN are of type 1, 



N N fN 

wg ) 2 (k 1 ,k 2 ,k 3 )= J ds 3 J ds 2 J ds, (,' k n ~ ■ 

fN fN 

(A2) 

We may also write Q (n> as sum of ordered integrals involv- 
ing different permutations of the wavevector arguments, 
of the form given in Eq. (|155[) , if it is understood that: i) 
The sum over permutations must be restricted to permu- 
tations that do not change any monomer types, but only 
that exchange the identities of monomers of the same 
type, and ii) The ordered integral <D I -™' ) i is defined by 
integrating over a subspace defined by the requirement 
that si < s 2 < ■ ■ ■ < s n and that the integral with re- 
spect to monomer index Sj may not extend beyond the 
block containing monomers of type ij. For example, 



,-,( 3 ) 



where 



ki,k 2 ,k 3 ) = ^2 ) 2 (k 1 ,k2,k 3 ) +w^ 2 (k 1) k 3 , k 2 ) 



(3) 



(A3) 



N s 3 fN 

u^ 22 (ki,k 2 ,k 3 ) = J ds 3 J ds 2 J d 



si ( e 



fN 



fN 



(A4) 

The value of an ordered integral cj^ i can be evalu- 
ated using a slight generalization of Eq. (|162p , in which 
the limits of integration in the ordered integral are in- 
terpreted in this manner, and in which the integrand is 
generalized to allow for the existence of different statis- 
tical segment lengths in different blocks. 
We also define functions 



-00 



(ki,---,k„) 



N -n Cj (n) 
(n) 



(k lr --,k n ) (A5) 



u^>^(k u ---,k n ) = N-^>_ ln {k ir --,k n ) (A6) 

Explicit expressions for these quantities may be writ- 
ten as functions of the 2n dimensionless wavevectors 
K ai = \c a biy/N/6, and of f± and / 2 . Because the ex- 
istence of two statistical segment lengths b\ and 6 2 pro- 
vides two ways to non-dimensionalize each wavevector, 
we will write u>ffi i and i for diblock copoly- 

mers as functions of ki , . . . , k„ , rather than as functions 
of dimensionless wavevectors. 



2. Three Point Correlations 

a. Homopolymer 

The dimensionless ordered integral u/ 3 ^ for a ho- 
mopolymer is given exactly, for arbitrary wavevector ar- 
guments, by 



lg(Kf)-g(K$) 



(A7) 



Each of the ordered integrals required in (|166[) may be 
evaluated using this general result. To obtain an asymp- 
totic expansion of each of the resulting integrals, we drop 
all terms that are proportional to e~Q (which are not 
analytic functions of 1/Q in the limit 1/Q — > 0), and ex- 
pand the remaining terms in powers of 1/Q- To 0(Q~ 4 ), 
this yields 



£(K,Q_,-Q+) 



1 

2Q2 



aK 
W 



w(Q_,K,-Q+) ~ -j 



2 
1 



K 2 (3 + 4a 2 )g(K 2 ) 



1 

(A8) 



An expression for u)(K,— Q + ,Q_) may be obtained by 
taking Q — > — Q, and a — > —a, in the expression for 
a>(K, Q,_ , — Q+). Terms linear in a cancel upon adding 
permutations, giving Eqs. (1173)1 and (|174|) . 



b. Diblock Copolymers 



For a diblock copolymer, we need the four functions 
~,(3)n, ~ ~ ^ ,-,(3), 

Hi 



^( k .q-^-q+) J ^j(k,q_,-q+), w$V(k,q_,-q+), 
and (k, q_, — q+) for i ^ j. The function 



(k, q_, — q + ) may be obtained by taking q — * q in 



our result for cD?- (k, q_, — q + ). The function Cb^/ for a 
diblock copolymer with an i block of length fi N is simply 
equal to the corresponding function for a homopolymer 
of type i and length faN. 

To calculate the remaining two functions, in which one 
of the species indices is different from the other two, we 
may start from the general result 



~,(3) 



£$(k ,k 6 ,k c ) = ei (K^)- 



Q-3 

K 2 - 



K 2 - 

a-3 



(A9) 



where K 2 ai = \k a \ 2 b 2 N/6 and and K% = |k c | 2 6fiV/6. To 
calculate £>f5, we note the symmetry w^-(k c , kb, k a ) = 

Or- - (k a ,kfc,k c ). Evaluating and expanding the required 
integrals yields 



)(?) 



(k,q_,-q + ) ~ BiiKDe^K 2 ) 



aK 



Q 2 i 



+ [e 3 (K 2 )K 2 (3 + 4a 2 )-A 



^(K 2 ) 
(A10) 



Expressions for (k, — q + , q_) and w ■ (k, — q+, q_) 
may be obtained by taking q — > — q in the above 
to integrals. The quantities Wjy-(q_, k, — q_|_) and 
w^. ■(— q + , k, q_) are C(Q -6 ), and so may be neglected. 
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Four Point Correlations 



a. Homopolymers 



The ordered integrals required to calculate the four 
point function ^ (K, -K, Q, -Q) may be expressed as 



VT'(K,±Q,tQ,-K) = J d 4 se- K2 ^H ± (s 32 )e 

o 

1 

^ (4) (K,±Q,-K, T Q) - J d 4 le- K2r ^H ± (s 32 )e-Q 2 ^ 

o 

l 

^ (4) (±Q,K,-K, T Q) - J dHe-Q 2 ^H ± (s 32 )e-Q 2 ^ 



where 



J ff±(s 32 ) = e- Q ±^ -e 



(All) 



and Q± = Q ± K. The dominant behavior of each of 
these intgrals can be obtained by noting that S32 is con- 
fined to very small values by the exponential factors in 
H(§32), and that we may thus approximate the integral 
with respect to S3 over the domain §2 < S3 < §4 to a first 
approximation by an integral over §2 < S3 < 00. This 
approximation yields a common factor 



00 

Jds 32 H ± ( S3 2) = ^T-q^ 



in the integrals 

V) (4) (k,±q,tQ,-k) 



^ 4 >(k,±q,-k,tQ) - i^g(K 2 



± ^- + K 2 

T2Ka 4(Ka) 2 
~Q^ + Q 4 



1 dg{K 2 
1 

2Q 



(A12) 



ds H±(s) 

/ dsH±(s) 
Jo 

1 r°° 

^ 4 )(±Q,K,-K, T Q) * 7^ / dsH ± (s) . (A13) 

v Jo 



Upon adding the four permutations in which ±K are the 
first and last arguments, and using the expansion given in 
Eq. (|Al2j) for the remaining integral, terms that are pro- 
portional to a and (D(Q~ 3 ) cancel. This leaves a leading 
order contribution of 0(Q~ 4 ). All other permutations 
lead to contributions of (D(Q~ e ) or higher, which do not 
lead to UV divergent contributions to 8uj^ in the one- 
loop approximation. 



$g,(k, _ k) q , _ q ) = N-*4%(k, -k, q, -q). As for ho- 

mopolymers, we may express fyjli (k, — k, q, — q) as a sum 
of ordered integrals. To calculate an ordered integral 

ip. . (ki, k.2, k3, k4) for a diblock, we require that the 
— 11 ■■■14 

integrals over each monomer index Si be constrained to 
the block specified by the corresponding monomer type 
index, as well as si < s 2 < s 3 < S4, and thus that 
ii < *2 < *3 < H- As for homopolymers, we find that 

-(4) 

tp^ kl (k 7 — k, q, — q) = for ordered integrals in which ±k 
are the first two or the last two wavevector arguments. 

As found for homopolymers, we find that the the dom- 
inant contributions to ip^li are 0{Q~ 4 ), and arise from 
ordered integrals in which ±k are the first and last ar- 
guments of For diblock copolymers, the (D(Q~ 4 ) 
contributions are obtained only from ordered integrals of 
the form V' ifcfc . (±k, ±q, =Fq, ±k), where ±q are associ- 
ated with monomers in the same block. The sum of the 
four ordered integrals of the form V^ 4 "^ (±k, ±q, =Fq, ±k) 

yield a contribution to V'niiXk, q, — q, — k) equal to that 
obtained from for homopolymer of length f\N . The re- 
sulting contribution to ^11x1 (k, q, q, — k) is 



,7,(4) _ WK 2 a 
Vim — ' 



8K 2 a 2 dgn 
Q 4 d(K 2 ) 



(A14) 



The dominant contributions to ^1112 (k, q, — q, — k) arise 
from the ordered integrals 7/> m2 (±k, ±q, =Fq, — k), which 
yield 







n 



where 



SK 2 a 2 8g 12 
Qt d(K 2 ) ' 



h - fi and s[ 



(A15) 



Cor- 



responding approximations for V>2222 O^k, ^li "FQj — k) 
7/^222 (±k, ±q, =Fq, —k) can be obtained by analogy, by 
switching the labelling of blocks 1 and 2. The dominant 
contributions are thus all of the form 



b. Diblock Copolymers 

Consider a diblock in which block 1 extends from 
< s < fxN and block 2 is f x N < s < N, and let 



,7.(4) „ % K W %J 
Vikkj- Q , ^ 



(A16) 
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